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Abstract 

At low Q^, charm production in deep-inelastic scattering is adequately described by assuming 
generation in electroweak boson-light parton scattering (dominantly boson-gluon fusion) which 
naturally incorporates the correct threshold behaviour. At high this description is inadequate, 
since it does not sum logs in Q"^ /rn^, and is replaced by the treatment of the charm quark as a light 
parton. We show how the problem of going from one description to the other can be solved in a 
satisfactory manner to all orders. The key ingredient is the constraint of matching the evolution of 
the physical structure function F2 order by order in as{Q'^) in addition to the matching of the value 
of F2 itself. This leads to new expressions for the coefficient functions associated with the charm 
parton which are unique in incorporating both the correct threshold and asymptotic behaviours at 
each order in perturbation theory. The use of these improved coefficients lead to an improvement 
in global fits and an excellent description of the observed F2^charm- 



RAL-TR-97-049 
September 1997 



1. Introduction. 

The factorization theory in QCD [1] has existed for many years, and has been one of the 
triumphs of quantum field theory. However, in its original form it does not take account of a 
number of possible complications, i.e. it exists only for massless particles, and its ordering does not 
take account of possible enhancements at high orders in at small x. Until the past few years 
both of these complications were not of any real phenomenological importance. The lowest values 
of X probed were large enough that it was unimportant to consider small x enhancement. Also, 
the up, down and strange quarks were considered light enough to be treated as massless whenever 
one was within the realm of perturbative QCD. Furthermore, there was little data on the charm 
contribution to the structure function and it was generally such a small component of the total 
structure function that it could be treated very approximately. 

Both of the above complications have recently become a great deal more important due to the 
advent of HERA. This now probes structure functions at far lower values of x than any previous 
experiments, as low as x ~ 10~^ [2] [3], and the treatment of structure functions should really take 
proper account of any small x complications. Also, the small x structure functions now have a 
contribution due to the charm structure function which is far from insignificant, i.e. it can be 
more than 20% of the total structure function and, moreover, in the past couple of years direct 
measurement of the charm structure function has also become possible [4] [5]. This has made it 
essential to treat the contribution to the structure function due to massive quarks in a correct 
manner. 

In this paper we propose a new method for the treatment of heavy quarks in structure functions. 
We begin by describing the features a correct treatment must exhibit at both high and low Q^, 
and the techniques used in either of these limits. We then give a discussion of the correct way 
to take account of heavy quarks in a well-ordered manner over the full range of Q^, showing how 
this relates to present treatments, and in particular demonstrating that one may choose to evolve 
the partons according to the massless evolution equations. However, we shall see that the detailed 
construction of the coefficient functions required is extremely difficult if not impossible. Therefore, 
we provide instead a prescription for calculating structure functions including heavy quark effects 
which is somewhat simpler than the strictly correct treatment, and which is directly analogous to 
the normal manner in which one calculates order by order for massless partons, but which is in 
practice essentially identical to the strictly correct treatment. Finally, we will present the results 
of a comparison of our method to data: both that for full structure functions and for the charm 
component to the structure functions. These comparisons turn out to be very good. We also 
make predictions for the charm component of the longitudinal structure function. Complications 
due to leading ln(l/x) terms at all orders in are ignored, and while a correct treatment of 
structure functions should of course deal with this problem, we feel that this would overcomplicate 
our presentation, and besides we wish to compare directly with normal NLO in approaches. A 
paper which takes account of both small x complications and massive partons is in preparation. 
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2. Structure Functions With Massive Quarks. 

We consider the case of n^- massless quarks and one heavy quark. One of the simplest ways 
to deal with heavy flavour production in deep-inelastic scattering is to treat the mass of the heavy 
quark, M, as a hard scale [6]. In this case the Uf light quarks are always treated as partons, but all 
other quarks are never treated as partons at any scale: the cross-section for production of heavy 
quarks is expressed entirely in terms of coefficient functions depending on the heavy quark mass 
convoluted with parton distributions which only depend on light partons, i.e. 

M') = CfJiCi'li^^ QVM^) ® ifi') + 0(AVm2), (2.1) 

a 

where (Jj(x, Q'^,M^) is the cross-section for scattering off a particular quark, either heavy or light, 
and a runs over the light partons, i.e. the gluon and the n/ light quarks. This approach is very 
well-defined in theoretical terms, essentially being a simple generalization of the usual factorization 
theorem, with (2.1) being valid to all orders up to the higher twist corrections of 0(A^/M^). 

This approach is adopted by a number of groups [7] , and is usually known as the fixed flavour 
number scheme (FFNS). It is normally used in the particular renormalization scheme where all 
diagrams with no heavy quark lines are renormalized in the MS scheme, while those with heavy 
quarks are renormalized at zero momentum. This is particularly convenient because the effect of the 
heavy particle decouples from the light sector, in particular the coupling is the 3 flavour MS coupling 
and the hght parton distributions evolve as in the 3 flavour MS scheme. The Cf^^{z, Q'^ / fi'^ , Q'^ / M^) 
have all been calculated to 0{as ) [6] [8] [9] and ©(a^) [10] in this scheme, though analytic expressions 
only exist at 0{as). 

In principle this approach is a very good way to calculate the effects of heavy quarks in structure 
functions. At each order it incorporates the kinematical threshold in the light parton-photon centre 
of mass energy = Q^{z~^ — ^) > 4M^ in a smooth manner (which then guarantees the same 
smooth threshold in the invariant mass of the hadronic remnant W^, up to proton mass corrections) 
and the coefficient functions are calculated order by order in precisely the same manner as the light 
particle coefficient functions (though the actual calculations are rather more difficult). However, it 
does have one major shortcoming. As one calculates to higher orders in ag one encounters higher 
powers of ln(Q^/M^) and ln(/x^/M^). Letting /x^ = Q^, and thus eliminating all logs in Q'^/jj,'^, 
then for — oo the coefficients at mt/i-order in (as(Q^)/27r)'" have the series expansion 

Cf/^"'{z,Q'/M') = Y,Uz)\n-{QVM^). (2.2) 

ra=0 

Thus, working order-by-order in ag in this approach one is failing to take account of these large 
logs. This is not only a practical concern in the sense that these large logs in Q^/M^ at higher 
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orders in Ug can potentially be phenomenologically important^, but is also a theoretical concern 
insofar as at each order in the leading power of ln((5^/M^) is the leading in ag part of the overall 
coefficient function with this ln(Q^/M^) behaviour, and is really part of the leading order expression 
for the structure function as a whole. The same reasoning applies for the next-to-leading power 
of ln(Q^/M^) etc.. This is similar in principle to the problem of increasing powers of ln(l/x) with 
increasing powers of . It is more difficult in one sense, in so much that in the expressions for the 
coefficient functions these large logs in (Q^/M^) are hidden within very complicated expressions. 
However, it is far simpler in the particular limit » because we know exactly how to sum 
the logarithms in Q^, i.e. we solve the renormalization group equation for fixed order in a^. 

Thus, in order to sum these large logs in Q^/M^ it is convenient to consider the heavy quark 
to be a parton and for its distribution function to satisfy the renormalization group (DGLAP) 
equations as do the light partons. An extremely simple approach which incorporates this idea is 
the zero-mass variable flavour number scheme (ZM-VFNS). This treats the massive parton as being 
infinitely massive below some threshold in /j,^, and totally massless above the threshold, i.e. all 
coefficient functions coupling directly to the charm quark turn on at the threshold, the evolution of 
the charm quark begins at this threshold, and the number of flavours in the coefficient functions, 
anomalous dimensions and the running coupling constant increases by one to n/ + 1 discontinuously 
at the threshold. Despite the simplicity of the approach this procedure must in principle be done 
with care if the correct results are to be obtained in the asymptotic limits [11] (see below for 
details). In particular, the decoupling theorem tells one how the coupling constant must change in 
order to get the correct results well below threshold. Also, the parton distributions just above the 
chosen threshold must be carefully defined in terms of those below threshold in order to guarantee 
that the correct result is obtained as — oo. In practice at low orders the situation is relatively 
simple, e.g., if the threshold is chosen to be precisely /x^ = M^, then at next to leading order, the 
light parton distributions are continuous across the threshold (in MS scheme) and the evolution of 
the charm parton distribution begins from zero. At higher orders the parton distributions must 
change discontinuously across the threshold and in particular the charm evolution must begin from 
a non-zero value. 

For many years the above approach was that most commonly used in global fits. The CTEQ 
collaboration used the approach at next to leading order, as explained above [12], while the MRS 
collaboration motivated their choice of threshold by phenomenological considerations rather than 
the strict theoretical ones [13], but in practice this resulted in a very similar choice of threshold (i.e. 
2.7GeV^ for MRS compared to 2.56GeV^ for CTEQ). While the charm contribution to the structure 
functions near the region of threshold was not too important this simple treatment was perfectly 

^ They are not important for « because the large logs are killed by factors coming from the 
kinematical threshold. 
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adequate. However, it is clear from its construction that it will not suffice as a good description 
of charm production in the region of the charm threshold. In particular charm production has a 
sharp threshold at a chosen rather than a smooth threshold in W"^. 

Hence, some approach which extrapolates smoothly from the FFNS at low to the ZM- 
VFNS at high is required in order to produce a good description of the effect of heavy quarks 
on structure functions over the whole range of Q^. Let us discuss how this may be achieved. In 
order to do this we first put the ZM-VFNS on a more solid theoretical footing. If we regard the 
quark mass as being a soft scale then the factorization theorem tells us that 

aM^M') = Y,Cl^+\Qy^^')^f:^+\^^^My^^') + 0{M'/^^% (2.3) 

b 

where b runs over the light partons and the massive quark. We are able to remove the large 
logs in Q^/M^ from the coefficient functions, and hence obtain the normal massless coefficient 
functions, and absorb them into the definition of the parton distributions at the expense of having 
potential "higher twist" corrections of 0{M^ //j?). While the parton distributions depend on M^, 
if the operators defining the partons undergo ultraviolet operator regularization in the MS scheme 
then their evolution depends only on the anomalous dimensions obtained from this ultraviolet 
regularization. These are independent of the mass of the heavy parton, and the evolution is as if 
for n/ + 1 massless quarks in the MS scheme. Hence we have the formal definition of the ZM-VFNS 
which will become exact for » M^. 

However, we have one more degree of freedom in (2.3) than in (2.1), i.e. we have the heavy 
parton distribution to parameterize at some arbitrary starting scale for evolution, and also no 
apparent reference to the mass-scale in the definition of the parton distributions. This is not 
in fact true since it can be shown that 

/;^+\.,/x^/.VM2) = A^'^(/xVm2)®/:^(/x2), (2.4) 

where the operator matrix elements A'"^{z, ^'^ /M^) contain logs in {/I'^/M^), and are calculable 
order by order in perturbation theory [14] [15]. Hence, the partons in the ZM-VFNS can in fact be 
generated from those in the FFNS at all /x^ by using the leading logarithmic expressions for the 
operator matrix elements and the expression (2.4), rather than using the four-flavour evolution 
equations at all. Indeed, if the starting scale is chosen as /i^ ^ then strictly speaking all the 
leading logs in (jl'^/M^) should be included in the matching condition which is just as complicated 
as using (2.4) at all scales. However, if the scale at which evolution begins is precisely /x^ = M^, 
then the matching condition for the partons in the two schemes is a power series in with no 
logs. Therefore, it simplest to use (2.4) only to define the order-by-order parton distributions at 
the starting scale, and then to calculate the parton distributions at other scales by evolving using 
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Uf + 1 massless flavours. This procedure guarantees the correctness of the ZM-VFNS calculation 
in the limit Q'^ » M'^. 

By comparing the expressions (2.3) and (2.1) at » M^, and using the relationship (2.4), 
one can calculate the FFNS coefficient functions, up to 0{M^/Q'^) corrections in terms of the 
massless MS coefficient functions for n/ + 1 flavours [14] [15], i.e. 

Cf/{z, Q^fx^Qyu') = Cr/+^(QV/.2) ^ A^'^ifiyu') + ^(MV/x^). (2.5) 

The detailed expressions of this form can be found in [15], where they are used to calculate the 
— > oo limit of the heavy quark coefficient functions in terms of the known light quark coefficient 
functions and calculated operator matrix elements. These authors then define f^^^^p as the 
structure function obtained from these asymptotic expressions for the coefficient functions and the 
parton distribution in (2.1). They then, through purely phenomenological motivation, define a 
variable flavour number scheme [15] [16] by the formal definition 

jpVFNS _ pZM-VFNS _ pASYMP _|_ pFFNS 

This then extrapolates smoothly from one limit to the other, being guaranteed to reduce to the 
correct limit order by order in as{Q'^) at high Q^, though only approximately to F^^^^ order by 
order at low Q^. 



(2.6) 



3. A Complete Treatment of Charm Mass Corrections. 

Although we agree with (2.5) and hence with the results at high regarding coefficient 
functions in [14] [15], we believe one may be more ambitious. Rather than simply accepting the 
uncertainty of 0(M^//x^) in (2.5) we can be more systematic and demand that there is a scheme 
which uses the definition of the parton distributions in (2.3) and (2.4) but which is correct up to 
0(A^/M^). Inserting (2.4) into (2.3) and subtracting from (2.1) it is clear that the difference is 

<(mV/^') ® = [cf/(gV/^^QVM^) - cl^^'iQVt^') ® ^''"(mVm^)] ® /^(m^), (3.1) 

where c"(z, M^/^u^) is 0{M'^ / ji^). Also making use of (2.4) this difference can be written as 

<(MV/x2) ® (A^«(MV/x2))-i /;^+^(/x2,/xVm2), (3.2) 

and so c1{M'^ /ij?) (g) (^^"(M^///^))~^ plays the role of the correction to the massless coefficient 
functions which accounts for the 0{M'^ / jj?) corrections. Defining a corrected coefficient function 
as 

CY{z,Q'/i^\QVm') = Cr/+\z,QVM') +<(MVa^^) (^^«(m7/x^))-i 

= (QVM^QVM^) ® {A'-{fiVM'))-\ ^^'^^ 
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then we have the factorization theorem 

a,ix, Q\M^) = ^ (QV/^^ MV/x^) ® K'^\^^^My^^') + 0(AVm2). (3.4) 

a 

Thus, (3.4) gives us a method for defining the structure function including a heavy quark up to 
errors of 0{h? [M"^) but where all partons evolve according to the massless evolution equations. It 
seems sensible that the best way to proceed for the calculation of structure functions in the presence 
of a heavy quark is to use the FFNS up to some scale of O(M^) and then switch to the scheme 
defined by (3.4) above this scale. (Changes of renormalization scheme across threshold applying 
to situations of this general type were first proposed in [17].) We shall call this a variable flavour 
number scheme (VFNS). Moreover, we believe that it is sensible to choose the renormalization 
and factorization scale /j,^ = in both schemes, for all scales and for both light and heavy 
quark structure functions.^ This very simple choice automatically avoids having different scales 
for different components of the complete structure function, and means that all mass effects are 
contained entirely within the coefficient functions. It also agrees with the normal asymptotic choice 
of /i^ = and removes all problems of logs of Q^//^^ (the solution of the evolution equations 
summing such terms) and /i^/M^, and we are left just with the problems of \n{Q'^ /M^)? This 
choice is expressed explicitly in all our equations from now on. Finally, as already mentioned, if 
we choose the transition scale as precisely fi^ = Q"^ = then all the logs in Q"^ /M"^ disappear, 
and the matching conditions between the partons in the two schemes in (2.4) are a simple power 
series in ag^M"^). Thus, performing the matching at M^, and solving order by order in OigiQ^) as 
in the strictly massless case we are guaranteed to sum the logs in Q"^ /M"^ correctly at zeroth order 
in M'^ IQ'^ . Combining with the mass corrected coefficient functions to the appropriate order, we 
should then get the mass corrected structure functions correctly order by order. Unfortunately, the 
procedure is not quite as simple as this. 

We see that the defining expression for {z, /M"^) is in fact of exactly the same form as 
(2.5), except that it is now exact at all Q^, rather than having corrections of 0{M^ /Q^), and that 
this time it is the n/ + 1 flavour coefficient functions which are the unknowns to be solved in terms of 
the FFNS coefficient functions and the operator matrix elements, rather than the asymptotic form 
of the FFNS coefficient functions. However, this leaves us with an ambiguity. The index a runs over 
the gluon and the light quarks while b also includes the heavy quark. Hence, while the asymptotic 
FFNS coefficient functions in (2.5) were defined uniquely in terms of the light n/ + 1 coefficient 

^ Of course, if we reach sufficiently low then we must introduce some finite renormalization scale in 
order to have a finite expression for heavy quark photoproduction. Since we only consider > IGeV^ we 
do not consider this problem in this paper. 

^ In the asymptotic expressions for the FFNS coefRcient functions in [14], this choice leads to significant 
simplification. 
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functions, solving (3.3) for the Cy^^{z, Q^/M'^) in terms of the FFNS coefficient functions does not 
lead to a unique solution. 

In order to demonstrate this let us write out our equations for the VFNS in full. For the case 
where the photon couples directly to the heavy quark, H, we have two equations: 

C'h7 = ^Lh ® CY/ + n,Al^^, ® C-f- + Ai, ® [cl^' + CY/'] (3-5) 

and 

^FF,5 _ .PS „ \f-,VF,NS , ^VF,P5-| , \ aNS , „ , aPS 1 ^ ^VF,PS , aS ^ ^VF,S /o (y^. 

^Hq - ^Hq ® [^HH + ^HH J + [^qq,H + nfA^^ jjl ® ^ Hq + ^gq,H ® ^Hg ' l^.bj 

where S, NS and PS refer to the flavour singlet, non-singlct and purc-singict (singlet minus non- 
singlet) respectively. In the case where the photon couples directly to a light quark we have three 
equations. Denoting the masslcss MS coefficient functions with Uf light flavours by Cia{nf) and the 
contributions to the light flavour coefficient functions in the FFNS due to heavy quark generation 
by C[g^ we have 

riNS I riPF^NS _ aNS ^ r'yF,NS (n ^\ 

^qq "T —^qq,H'^^qq ) ' j 

,^FF,S _ aS ^^VF,NS,aS i^ryF,S aS i^ryF,PS,AS ^(^VF,PS /o ON 

^qg + - '^qg,H ^^^qq + ^gg,H ® qg + ^f\g,H ®^qq + ^Hg ® qH l^.Sj 

and 

^qq +^qq - [^qq,H +n f A^^ jjI ®C ' +A^^(g)C^^ +Aqq^H®(^qq +Agq^H®(^qg ■ 

(3.9) 

These arc very similar to the equations (2.31)-(2.35) in [15] and, as in those equations, we have 
implicitly divided all pure singlet quantities coupling to quarks and all singlet quantities coupling to 
gluons by Uf. Also, as in these previous equations, it is implicit that all quantities on the left-hand 
side arc expanded in the n/— flavour MS coupling constant while those on the right-hand side are 
expanded in terms of the (n/ + l)-flavour MS coupling. The relationship between the two couplings 
was calculated in [18] and corrected in [19]. It is 

a,,„,+i(Q2) =as,nf{Q'')+al^,{Q^)^Tf\n{Q''/M^) 

1 r 1 1 1 9 1 

+ [g^l \n\Q^/M^) + ^{^CATf + ACfTf) HQ^M^) + ^TfCf - -T;C^] + • • • , 

(3.10) 

where the coefficient of the leading log at each order in as,nj(Q^) is the same in all schemes, but 
other coefficients depend on details of renormalization, in particular whether the mass M is the 
fixed or running mass. The particular choice above corresponds to a fixed heavy quark mass at 
NLO. 

The difference between our expressions for the coefficient functions and those in [15] is that 
the coefficient functions on the right-hand side are the VFNS coefficient functions. Not only does 
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this mean that the equations are meant to hold including terms of 0{M^ /Q"^), and that we solve 
for the coefficient functions on the right-hand side, but also that there is a difference between 
the coefficient functions which couple to the heavy quark distribution and those coupling to the 
light quark distributions. For example, while C^^^^ and C^^'^^ must be identical in the limit 
Q"^ ^ CO they certainly do not have to be identical at moderate Q^, and physical intuition suggests 
they should not be. This means that unlike [15] wc do not have five equations for five unknowns, 
but we have five equations for eight unknowns. In order to reduce to the correct ZM-VFNS at 
very high we must choose definitions for the mass-corrected coefficient functions which reduce 
to the n/ + 1 light parton coefficient functions as — >^ cxd, but this constraint still leaves a great 
deal of freedom. 

As an example let us consider what is in practice the most important case, the equation for the 
boson-gluon fusion coefficient function for the heavy quark structure function F2,i/(x, Q^), (3.5). 
The expansion of Cl^^ begins at 0{a,{Q'^)) as does C^^^^ and ^1^, while A^^^ jj and C^^^^^ 
begin at zeroth order. Using the known expressions for the operator matrix elements we obtain the 
lowest order equation relating the FFNS coefficient functions and the VFNS coefficient functions 

<Hf '(^, QVM^) = cZf\^, Q'/M^) + MQ^M') + c.,)P°, <5'^^'°(QVM^), (3.11) 

where Pqg{z) is the lowest order splitting function, and Crs is renormalization scheme dependent, 
but Crs = in MS scheme. Hence, we have freedom in how we choose our zeroth order heavy 
quark non-singlet coefficient function, and this then determines our first-order mass-corrected 
gluon coefficient function. More generally, we have freedom in how we define each of the three 
coefficient functions coupling to the heavy quark C^^^^ , C^^^^ and C^^ '^^ each order in 
perturbation theory, being constrained only by the requirement that they are of the correct form 
as oo. 

Of course, there cannot truly be an ambiguity in the order-by-order definition of the structure 
functions. In order to illustrate this consider the structure function F2{x,Q'^). We also come back 
to the point concerning renormalization scheme dependence. In order to maintain renormalization 
scheme consistency we must be very careful about the way in which we order the expressions. 
Doing this correctly does not remove the ambiguity in our definitions of the coefficient functions, 
but it does render this ambiguity physically meaningless, even order by order. Let us consider 
specifically the heavy quark contribution to the structure function F2{x, Q^) in the general VFNS. 
In fact we will discuss its In Q^-derivative since it is the evolution of F2^h{x,Q^) which is a more 
natural quantity. Taking the InQ"^ derivative of F2^h{x, Q^) and keeping all terms up to 0{as{Q'^)) 
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multiplying the VFNS parton distributions we obtain 



dln(Q2) dln(Q2) 

^VF,NS,0 



{H{Q') + H{Q')), 



+ 



27r 



2,HH 



+ 



27r 



+ 



{Q''/M^) 

dln(Q2) 
dln(Q2) 



(3.12) 



Asymptotically the second and third term in this expression reduce to the required form for the 
leading order expression in the ZM-VFNS. All other terms fall off to zero in this limit so we are 
guaranteed the correct asymptotic expression using this prescription. However, at low the 
exact form of the expression is highly sensitive to our particular choice of coefficient functions. 
This clearly means that we do not have a truly well-ordered solution and this is because the true 
ordering of the coefficient functions C2fji,{z,Q'^/A!^) is not as simple as just order by order in 
as,nf+i{Q'^) due to their dependence on the quark mass. Indeed, their ordering is the crux of the 
problem, and we will explore this below. 

In order to examine the true ordering of our expression we will express it in terms of unambigu- 
ously defined quantities, and also in terms of those where the ordering is relatively straightforward. 
Hence we will we will express it in terms of the FFNS parton distributions, the mass-dependent 
coefficient functions CYh},{z,Q^ /M^), the operator matrix elements and the coupling ag^ufiQ^)- 
The FFNS parton distributions are correctly ordered simply by solving their evolution equations 
to a given order. The operator matrix elements are ordered according to the power of as^uf+iiQ"^) 
minus the power of InlQ"^ /M"^), i.e. the leading-order term is of the form 



n=l ^ ' 



)an{z). 



(3.13) 



The n/ + 1-flavour coupling constant is defined in terms of the n /-flavour coupling in an analogous 
manner, i.e. the leading order relation is 

a,,r.,+im = a,,^,+Y.<n){Q\-^) \n-{Q^/M% (3.14) 



First using the expression for C^^'giz, Q^/M^) (3.11), but only keeping the leading order part 
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of the operator matrix element, i.e. leaving out the Crs, and substituting into (3.12) we obtain 



dln(Q2 
+ 

+ 



{H{Q') + H{Q'))o 



) dln(g2) 

Q:g,n/ + l (Q^) ^VF,NS,0/ 
as,n^ + l(Q^) 



27r 



dln(Q2) 



+ 



dln(Q2) 



9S 



27r 



dln(Q2) 



(3.15) 

We can then be more detailed by using the explicit expressions for {H(z, Q'^) + H(z, Q^))o and 
gQ^^^{z, Q^) in terms of the FFNS parton distributions, i.e. 



(^s,nf (Q 



ln(QVM2)n>5o'(Q') + ^«n,(Q')ln^(Q7M^)), (3.16) 



2\ i„2/^2 /nyr2\ 



27r 



{H{z,Q-') + Hiz,Q-'))o-^ 
and 

9o'^\z, Q') = g7 (z, Q') - ^^^^^ HQVm')97 {z, Q') + iQ') In^Q'/M^)), (3.17) 

and also the expression for the n/-flavour coupling, (3.14). Doing this and remembering that 

da 



dln(Q^) 



C>(M2/Q2) then we obtain 



ciln(Q2) 



+ 



27r 

<-'^s,n/(Q ) 

2^r 

r2 //o2> 



ciln(g2) 



ln(QVA^^ 



' dln(Q2) 



(3.18) 



ln(QVM^)<^'^^'°(QVM^) P°, P°, 5o ^ (Q7 

+ 0(M2/Q2) . 0(a2^^^(Q2) 1^2(^2/^2)^ 



Hence, as well as asymptotically reducing to the correct leading order expression, the prescription of 
keeping all terms up to 0{as^nf+iiQ'^)) which multiply the leading order VFNS parton distributions 
has resulted in a unique 0{aa^nfiQ'^)) expression for the derivative of the heavy quark coefficient 
function which also (and necessarily) has the correct threshold behaviour. However, it is clear that 
the C'(a2^^ (Q^) InQ^/m^) expression, while having the correct asymptotic limit, has behaviour 
for ~ M2 which is sensitive to our choice of coefficient functions. In particular, the behaviour 
of these terms will not generally respect the threshold in W^. It is clear that at higher orders 
in as,n/ (Q^) ln(Q^/M2) while we will obtain the correct asymptotic behaviour, our lowish Q'^ 
behaviour will be dependent on the choice of coefficient functions. 
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If we were to use the expression for the structure function itself, rather than its derivative, in 
the VFNS by combining the lowest order in a^^nf+iiQ^) coefHcient function with the lowest order 
VFNS parton distributions, i.e. 

F2,Hix, Q2) = CI^hh''\QVM^) C5 mQ^) + H{Q^))o, (3.19) 

then again we would be guaranteed the correct LO expression in the asymptotic limit. However, 
even the leading term in (Q^) ln((5^/M^) (when expressed in terms of the FFNS parton 
distributions and operator matrix elements) is now completely dependent on the choice of coefficient 
function, and there is no requirement to have the correct threshold behaviour at all. 

It should be no surprise that we have this problem. As mentioned earlier in the FFNS the 
coefficient functions to all orders contain renormalization-scheme-independent leading-order con- 
tributions. By working in the VFNS we have managed to extract the asymptotic form of this 
leading order contribution in a relatively simple manner. However, in order to have the full leading 
order expression for the structure functions in the VFNS in the threshold region, we need to extract 
all the information from the leading order contribution to the FFNS coefficient functions. In prin- 
ciple, the full LO FFNS expression should contain the leading parts of the coefficient functions at 
all orders in as^njiQ"^)-, and the LO VFNS should include coefficient functions constructed from the 
full LO FFNS coefficient functions and the full LO operator matrix elements. Absolutely correct 
matching between the FFNS and the VFNS at = leads to the absolutely correct renor- 
malization scheme consistent description of both of these schemes. Thus, in practice the strictly 
correct LO VFNS is no simpler than using the strictly correct LO FFNS coefficient functions. This 
is extremely difficult indeed, and in fact probably impossible, there being no clear unique way in 
which we subtract out the leading order, renormalization scheme invariant part of the ©(a" „^ (Q^)) 
FFNS coefficient function except in the asymptotic limit. Indeed, if we were to proceed further 
for our above example of '^^^\^\^q?^ ^ we would find that our definition of the LO contribution at 
^(•^s.n/ (Q^)) would rely on being able to extract an unambiguous LO, renormalization scheme 

independent part out of — ^'^rfin(Q^) ■ Though this is simple in the limit — >■ oo [14] [15], 

there does not seem to be any good prescription for arbitrary Q"^ . Therefore it appears as though 
the VFNS is only any advantage at all in so much that it gives a definition of the charm parton 
distribution. There does not seem to be any tractable way to produce a prescription for calculating 
heavy quark structure functions which both correctly sums the leading logarithms and which has 
absolutely correct, unique threshold behaviour. 
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4. A practical VFNS. 

Bearing in mind the difficulty, or indeed probable impossibility of producing the unambiguous 
well-ordered calculation of structure functions, it is our aim to produce a prescription for calculating 
heavy quark structure functions order by order in ag^Q"^) in such a way that we obtain relatively 
simple expressions yet maintain as much accuracy as possible over the whole range of Q^. Let us 
first consider the region of = and below. In this case if we work order by order in a^^nj (Q^) 
in the FFNS, i.e. define the rat^-order expression for the heavy quark structure function by 

F^^hM') = E E(^^^#^)"""^'<5r"^^^(MVQ^) ® /^^(Q^), n = 1 ^ oo, (4.1) 

we know that the strictly leading order terms we ignore are really an order of as,n/(Q^) down 
on those we keep, with no large ln((5^/M^) enhancement, for these values of Q^. Adopting this 
procedure, when working to 0(q;"„^(Q^)) we have an error of 0(q;"+^ (Q^)) compared to the (in 
principle) correct calculation, which is the same size as terms not yet calculated and the same size 
as the renormalization scheme uncertainty. This seems perfectly satisfactory for this region. 

Above = we want to order our calculation as in the massless case so that in the asymp- 
totic limit of >> we will obtain correctly ordered expressions. Therefore, we order the 
calculation by using up to 0(a"„j,+i(Q^)) coefficient functions when solving the evolution equa- 
tions using up to C'(a"^J_|_^((5^)) anomalous dimensions, as required by renormalization scheme 
consistency, e.g. the leading order expression is 

FU^,Q')=j:c2,^"(^yQ")^fo'^'(Q")^ (4-2) 

a 

the next-to-leading order expression is 

a 

+cr'°(MVQ^)®/r^+'(Q')], 



(4.3) 



etc.. We stress that this is not a choice, but a strict requirement of obtaining ordered asymptotic 
expressions for the structure function itself or its ln{Q'^) derivative. Of course, in this region 
of we now have the ambiguity in the definition of the coefficient functions. Thus, since we 
are not performing the strictly correct ordering we have to make a choice for these coefficient 
functions. We do this by defining them order by order in as,n/+i(Q^) using the equations (3.5)- 
(3.9) which guarantee correctness to all orders, and also by using the freedom to choose some 
coefficient functions, i.e. the three coefficient functions coupling to the heavy quarks, to bring us 
as close to the really correct calculation as possible. 

In perturbation theory it is not really the structure function at a particular value of for 
which we solve but the evolution at all in terms of the structure functions at some particular . 
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Bearing this in mind it seems sensible to constrain our coefficient functions by making tfie slope 
of the structure functions at a given order in as{Q^) to be continuous across the transition point. 
In order to examine this constraint, let us again consider the In Q^-derivative of F2^h{x,Q^)- 
Approaching the transition point from below, our prescription gives the lowest order as,nf{Q^) 
expression for the ln((5^)-derivative as 

c/ln(Q2) - 27r dHQ^) ^ >■ ^ > 

Just above = the ln((5^)-derivative of the LO expression in the VFNS is 

dln(Q2) ^ 

(4.5) 

where at the transition point the coupling at this order is continuous. Also we see that the artificial 
zeroth order term in (4.5) disappears at = (it is actually cancelled in the complete calculation 
as seen in (3.12)-(3.18)), and it is indeed possible to demand the continuity of the derivative 
across the transition point. Using the constraint and our simple prescription for constructing the 
structure function in the two regions we now have a unique form for the previously ambiguous 
C^^HH^'^{z, Q'^/M^). Using the fact that {H{z, Q^)+H{z, Q^))o = at = m^, we immediately 
obtain 

at = M^, and we define C2hh^'^{z, Q^/M"^) by demanding that it satisfy this relationship at all 
. As well as guaranteeing the continuity of the evolution of the structure function this definition 

also reduces to the correct form for » since in this limit — ^'"din{Q'^) ^ ^qg(^^ 

we shall see explicitly in §5). This means that the evolution will clearly reduce to the correct 
asymptotic form of a delta function in the limit Q'^ ^ oo. Above = terms are not exactly 
as prescribed by the absolutely correct procedure explained in the last section but they do explicitly 

maintain the correct threshold behavior since — '^'"d\n(Q'^) ^^^^ ^'-'^ ^ 4M^. At leading 

order we have in principle an error of 0{a^{Q^)) at the transition point due to the truncation of the 
FFNS expansion at 0{as,nf{Q^)) (where this error falls like (M'^/Q^) as we approach the correct 
asymptotic limit) and an error generated by the evolution which is zero at the transition point, and 
grows like as,nf+i{Q^) ln"((5^/M^), but falls like {M^/Q^) as we evolve up from this point. These 
errors are quite minimal, always being small compared to the quantity being calculated. Prom 
(3.11) we see that we have also completely defined C^^'j^'^'^(z, Q^/M^), i.e. in the MS scheme it 
is 

QVM^) = QVM^) - 1-(QVM^)^^^§|^^, (4.7) 
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though we have not yet made use of this coefficient function. However, we notice that each term 
in this coefficient function separately has the correct threshold behavior in W^. 

At leading order in this prescription the effect discussed above is the only real complication, 
i.e. the choice for C^^'^^'°{z,Q'^ /M'^) is the only one to make. Above the transition point the 
evolution equations for the partons are now in terms of n/ + 1 massless quarks, and the coupling 
constant becomes the MS coupling for n/ + 1 massless flavours. But all parton distributions and 
all other zeroth order coefficient functions are continuous across the transition. 

Of course, although we have determined the lowest order derivative of the coefficient functions 
on both sides of the boundary we must also discuss the value of the structure function itself at 
= M^. Using the zeroth order expression (4.6), the vanishing of the charm quark distribution at 
= leads to the charm structure function being zero there. Likewise the fact that at zeroth 
order in the FFNS the coefficient functions for charm production all vanish leads to the zeroth 
order value of F2^h{x, M^) being zero also. Thus, the two expressions are consistent. However, 
this is unsatisfactory for two reasons. Firstly, the leading order (order as{Q'^)) derivative of the 
charm structure function is non-zero both above and below = M^, provided x is low enough 
that we are above the threshold in W^. Hence, starting with a value of F2^Hix, M'^) = would 
lead to negative values for this structure function for < M^. Also, one would naturally expect 
the LO expression for a quantity to be a reasonable approximation to the quantity itself. The value 
of F2^h{x-,M'^) is not zero, and so the zeroth order expression is not a good representation of the 
true value. These problems come about because of a peculiarity of F2{x,Q'^) already discussed in 
[20]. In general its value at a given begins at zeroth order in (Qg), but the 0{as{Ql)) term 
is also really part of the leading order expression since it is rcnormalization-scheme independent. 
In contrast the derivative begins at 0{as{Q'^))., and all corrections are renormalization scheme- 
dependent and genuinely higher order. Thus, as argued in [20], the input should contain both the 
zeroth order term and the 0{as{Qf^)) term, but the latter should play no part in the evolution. 

Adopting this procedure we can now specify our leading order expressions for the charm struc- 
ture function as follows. Below Q"^ = we take the LO expression to be 

F2TM') = ^^^^<S(QVM^)0 5o^(Q^), (4.8) 

which is equal to the 0{as) value at = M^, and incorporates the LO evolution down from this 
scale (up to small corrections). Above = the LO expression is 

<^'"(^,Q') = Fl^'\x,M^) + (QVM^) ® {H{Q') + H{Q'))o, (4.9) 

which (up to the constant term) is of the standard form (4.2), and incorporates the correct LO 
evolution. In practice the constant term becomes almost insignificant as soon as > 4M^. Now 
we should consider the NLO expressions. 
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At next to leading order the situation is rather more comphcated because more terms come 
into play. We now define FFNS expressions by including terms up to order a^^nfiQ^) relative to 
the lowest order parton distributions. The NLO VFNS expression is defined as in (4.3) (up to a 
constant again). At this order the situation becomes more complicated because the pure singlet 
FFNS coefficient function becomes non-zero as does the contribution due to coefficient functions 
where the photon couples to a light quark but where heavy quarks are generated. Let us examine 
the NLO expressions for the derivative of the heavy quark structure function. First consider the 
C>(a^ „^ (Q^))) expression for the derivative of the heavy quark structure function in the FFNS. 
This is 



dF2,H{x,Q^ 
dln(g2) 



27r 

iFF.l, 



+ 



+ 



'2.Hg ^ 



dln(Q2) 



(4.10) 



a.n,(Q^)<g(QVM2) 



27r 



where {z,Q'^) is the singlet light quark distribution.^ In the VFNS the situation is even more 
complicated. Taking the derivative of the NLO expression, and ignoring those terms already in 
(4.5), we obtain 



dln(Q2) 



dln(Q2) 

«s,n/ + l(Q^ 



{H{Q^) + H{Q^)), 



+ 



27r 

+ cIZ{qVm') 



®{H{Q^)+H{Q^))q + 



dCl'/AQ' 



+ 



'2,HH 

~iVF,l 
'2,HH 



{H{Q^) + H{Q^))o) 

I {H{Q') + H{Q'))o + ^ g7^' {Q')) 



P 



r(Q7^ 



(4.11) 



^ We label Pgg(z) by the number of flavours because it is the only leading-order splitting function which 
depends on this number. The decrease of this splitting function above a threshold accounts for the fact 
that there is a new parton distribution, and guarantees overall conservation of momentum in the evolution. 
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These expressions are very difficult to compare in general. However, expressing the four flavour 
quantities in terms of the three flavour quantities the two are identical at NLO in MS scheme at 
= (the discontinuities in both the parton distributions and the coupling begin at NNLO). 
Thus, the heavy parton distributions {H{z,Q'^) + H{z,Q^))o and {H{z,Q^) + H{z,Q'^))i vanish 
at this point, and so do many other terms in (4.11). Prom the definition of (^^^-'^(-2, Q^) we can 
see that the term depending on gi{z, Q^) is the same in both expressions, and using the equation 
(4.7) we can see that '^'^d^°Q2')^ ^ = at = M^. Also in the combination -Pnff{z)+P9!p^ ^ f 
the flavour dependence cancels between the two terms, so this combination is the same in both 
expressions at = M^. 

Thus we have a great deal of simplification when comparing the two expressions at = . 
As in the LO case we can equate the terms coupling to the gluon in the two expressions, i.e. 



''''''fln(Q^^'^ = <^'1(QVM^) + CZ'UQVM^) ^ P^r^\ (4.12) 



at = M^, and this serves as a definition of the coefficient function C^^'^(2:, Q^/M^) at this Q^. 
However, unlike the LO case wc cannot define the coefficient function at all simply by extending 
this expression to all Q^. This is because it will not result in the correct asymptotic expression for 
C^Jiniz^Q'^ /W^), i.e. — '^'"iia(Q2) — - contains a hi{Q'^/M'^) term which must be cancelled. It is 
quite easy to find the generalization of (4.12), however. If one differentiates both sides of (3.5), and 
keeps those terms of ©(a^^^^ (Q^))) which survive as ^ oo (all terms of the form '^'''^'''^inlQ'^)^ ^ 
vanish in this limit since the VFNS coefficient functions tend to constants) then one obtains 

+ <^'^^'°(QVM^) ® ""^^ffg^f ^ + I HQVM-)C-Z{Q^/M-) « 

(4.13) 

where the last term comes about from the difference in the derivatives of the three and four flavour 
couplings. This expression guarantees the correct asymptotic expression for C^jJ'^(z, Q^/M^), 
while (4.12) guarantees the continuity of the NLO derivative of F2^h{x,Q'^) in the gluon sector, 
and hence the definition of C^^iz, Q^/M^) must satisfy (4.12) at = and (4.13) as oo. 
In fact, at = the two expressions are identical, i.e. 

dln(g2) ' 19 

and 



Pa.iz), (4.14) 
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and we have the very neat result that (4.13) is the generalization of (4.12) for all Q^, and 
C^Jilj^z^Q"^ /M^) is defined by (4.13). 

The above definition of C^Hui^i Q^/-^^)) when substituted into (3.5), determines the expres- 
sion for C^^'^(z, Q^/M^) which will be used at NNLO. However, we have now used up our single 
degree of freedom involved with the heavy quark structure function at NLO. Looking at the terms 
coupling to the singlet quark distribution in the two expressions (4.10) and (4.11) we find that the 
first contains 

cI'hUQVM^) ® + ^ (4-16) 

while the second contains 

Cr/AQ'/M') ^ + C[ZiQ'/M') (4.17) 

There is no degree of freedom in either of these equations, and no reason for them to be equal at 
= M^, and they are not. Indeed there was no further degree of freedom in the relationships 
(3.5) and (3.6) required of the heavy quark coefficient functions. Up to this order the only one 
available was for C2^^^(z, Q^/M^), and this has been determined by imposing the continuity of 
the evolution of the structure function in the gluon sector. Indeed, looking at (3.6) at 0{al{Q'^)) 
we see that we have already determined C^^'^^'^{z,Q^ /M'^), i.e. 

Cl^'r'\^, Q'/M') = Cl^n''\z, Q'/M^) - €^^{0" jM^) ^Jf (QV^^). (4.18) 

Using the framework wc have chosen to define the structure functions this discontinuity in the 
derivative of the heavy quark structure function in the singlet sector is unavoidable. There are 
simply not enough degrees of freedom to avoid it. In practice, since the evolution of the heavy 
quark structure function is driven very largely by the gluon, since this discontinuity begins only at 
NLO, and since (4.16) and (4.17) are not too different at = the effect is tiny. Of course, 
any discontinuity is only an artifact of the manner in which we are forced to do our fixed order 
calculations, and would disappear if we were to work all orders. In fact one can show that the 
discontinuity of the derivative in the singlet sector gets formally smaller as one works to higher 
orders. 

So now we have the definition of our NLO expressions for the heavy-quark structure function 
both above and below threshold. In the FFNS the definition is the simple extension of (4.8), being 
just 



= ^^^^(<S(QVM^)05o^(Q^) + <g(QVM^)05r(Q^)) 

(4.19) 



2\ 

2tt 
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which is equal to the 0{al) (i.e. NLO) value for the structure function at = and incorporates 
the NLO (i.e. 0{a^,^^{Q^)) evolution down from this scale (up to small corrections). The VFNS 
NLO expression is 

Q') ={^^) ' ^ 9o' (M^) + C7Si(l) So^ (M^)) 

+ CI^hUQVM^) ® {H{Q^) + H{Q^)), (4.20) 

which again, up to the constant term, which is the NLO input (the LO part of the input now being 
included automatically) , is of the standard form and incorporates the correct NLO evolution across 
the transition point. 

At this order we have to make some decision about how we treat the light quark sector. The 
lowest order contribution the heavy quark makes to a light sector FFNS coefficient function is for 
the nonsinglet coefficient function at 0{a^,^^{Q^)). Thus in the matching conditions between the 
FFNS coefficient functions and those in the VFNS in the light quark sector there are no mass- 
dependent corrections to the VFNS coefficient functions to 0{as,nf+i{Q^))- Hence, the evolution 
of the light quark coefficient functions above = A'P is exactly as in the massless Uf + 1 flavour 
case. Nevertheless, wc must decide on the form of the structure function at Q"^ = Af ^ and below 
this transition point. For the heavy quark structure function wc have been keeping heavy quark 
coefficient functions to one order higher in as{Q'^) in the FFNS than in the VFNS. This has been 
for the reason that the explicit ln(Q^)-dependence in the coefficient functions means that they 
contribute to the ln((5^)^derivative of the structure function at effectively one higher order in 
dsiQ^) than the VFNS coefficient functions, and also because the lack of the usual zeroth order 
coefficient function makes the ©(og^^^ (Q^)) coefficient function the LO one, the 0{a1 ^^{Q"^)) 
the NLO one, etc.. For the light structure functions there is a zeroth order coefficient function, 
so the second argument no longer holds. However, the former one still does, i.e. differentiating 
the expression for the light quark structure function below = and keeping terms of order 

'^s,n/(Q^) then — ^''^'^ riin(Q^) appears in the expression. This contribution accounts for the 

effect of the heavy quark to the evolution turning on as increases. For this reason we continue 
to keep the coefficient functions containing heavy quarks to one higher order than those with only 
light quarks even in the light sector. 

For the heavy quark structure function, because we had terms of higher order in q;s(Q^) below 
than above it, in order to impose continuity of the structure function at = we had to put 
a contribution to the VFNS expression which is constant, and one order in higher than the rest 
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of the expression (we also justified this from renormahzation scheme consistency). We now have to 
adopt a similar procedure for the light quark expressions. The NLO expression for the nonsinglet 
structure function for Q"^ < is 

(4.21) 

That for > is equal to 



+ ^f-d_imc-/-/+M ^ ;;^-.+i(Q2) ^ (^^y<-'--'^(l) /o"^-^(M^)). 

(4.22) 

In principle both sides should also contain a term oc a^((5o) for the genuinely light NLO input, 
where Qq is the scale at which the inputs are chosen. Such a term is always ignored, and would be 
very small. In practice all the 0{a1{Q'^)) terms in the above expression axe extremely small as well. 
The 0{a1{Q'^)) evolution derived from the above equations is not precisely continuous at = 

due to terms of inverse powers of Q"^ /M"^ present in — ^''''^ rfin(Q^) • '^'^^^ discontinuity will 

decrease as we go to higher orders, and these mass-dependent terms get absorbed by higher order 
mass-dependent VFNS coefficient functions. We note that leaving the 0(a^ „^(Q^)) term out of 
of (4.21) would also lead to a discontinuous evolution (actually more so) since the evolution would 
take account of tt,/ massless flavours below threshold but rz/ + 1 massless flavours above threshold. 

Finally at NLO the light quark pure singlet structure functions have no complications due 
to the heavy quarks at all. The flrst non-zero FFNS coefficient functions do not appear until 
0(Q;f „^((5^)), and so do not contribute to the evolution until NNLO. So at NLO we just use the 
n f massless flavour expressions below = and the rif + 1 massless flavour expressions above 
= M^. Continuity of both the structure function and its evolution are automatic. 

One could in principle work to progressively higher orders, but of course in practice the NNLO 
splitting functions and the NNLO FFNS coeflicient functions are all unknown at present. Neverthe- 
less, we outline the procedure to be adopted at all orders. For the heavy quark structure function 
there is essentially nothing new as we progress to higher orders. At nth nontrivial order we include 
all FFNS coefficient functions up to order as,nfiQ'^)^ ^^'^ VFNS coefficient functions up to 
order <-J+i(Q2). in the VFNS expression we always include the C'(a"(M^)) term which ensures 
continuity of the structure function. We determine C;^^'^~^(2;, Q^/M^) by demanding continuity 
of the derivative of the structure function at C(q!") in the gluon sector, and this determination 
predetermines ClH'g{z,Q^/M^) and C2^^'^"(z, QV^^) by using (3.5) and (3.6) to 0«_„^(Q2)). 
At C'(Q;s,nj,+i(Q^)) the coefficient function C^^'^(z, Q^/M^) becomes the sum of the nonsinglet 
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and pure singlet coefficient functions. Neither the condition (3.5) nor the continuity of the structure 
function and its derivatives determine these two contributions separately, so we are free to separate 
them as we wish, using the condition that each tends to the correct asymptotic limit. It would also 
be desirable to choose each so that they respect the kinematic threshold. 

For the light quark structure function the procedure at higher orders is also straightforward. At 
nth nontrivial order we include all pure light quark contributions to coefficient functions below = 

up to order ot^^^^iQ'^)', all mass dependent FFNS coefficient functions up to order Q:"„^((5^), 
and all VFNS coefficient functions up to order Q;"^J+i(<5^)- In the VFNS expression we always 
include the 0(q;"(M^)) term which ensures continuity of the structure function. Starting with the 
^(o^s,n/+i(Q^)) coefficient function we determine C^^~^(2:, Q^/M^) by demanding continuity 
of the derivative of the light structure function at 0{ag{Q^)) in the gluon sector, analogously 
to the heavy quark sector. With this one degree of freedom eliminated in this way all other 
VFNS coefficient functions are determined uniquely order by order in as by (3.7)-(3.9), i.e. this 
determination ofClg^~\z, Q^/M^) predetermines Cl^f{z, Q^/M^) and Cl^f^'^'iz, Q^/M^) by 
using (3.8) and (3.9) to 0«(Q2)). 

Thus, we have completely defined our prescription for calculating order by order for the struc- 
ture function F2{x,Q^). We can sum it up in the form of a table. This is shown in table 1. 
This method uniquely determines all VFNS coefficient functions, and while not leading to abso- 
lutely correctly ordered expressions it is a relatively simple prescription for obtaining order by order 
structure functions which are very similar to the strictly correctly ordered ones, which reduce to the 
correctly ordered expressions in the asymptotic limit and which order by order are consistent with 
all physical requirements. All prescriptions which obey (3.5)-(3.9) will be correct when summed to 
all orders, but some ways of choosing the heavy quark coefficient functions will clearly stay closer to 
the correct ordering than others. We believe that our prescription is the best available at present, 
and we see no easy way to improve upon it. 

We will demonstrate the results using our prescription in the next section, and see that indeed 
they do seem to work very well. However, first let us mention another currently available VFNS, the 
ACOT scheme [21] [22] [23]. Although there is currently no all-orders [24], or even NLO definition 
(for developments see [25], of the ACOT VFNS (which we will denote by ACOT) in print we 
believe that the definition of the coefficient functions in this scheme must be equivalent to that in 
(3.4), i.e. the VFNS coefficient functions are related to those in the FFNS by the equations (3.5)- 
(3.9). Indeed, at what they call LO, the ACOT coefficient functions satisfy (3.11). However, they 
determine the expression for C^^'h^'^ from the tree-level diagram for a massive quark scattering 
from a boson, and for a photon this gives 

C];,r'%,QVM^) = z6ixo-z)(l + '-^y fo= + (4.23) 
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where the modified argument of the delta— function follows from demanding the on— shell condition 
for the massive quark, and the remaining factor follows from the parton model for the longitudinal 
structure function, Fl = AM'^/Q^, which is added to the transverse component to give F2. Inserting 
into (3.11) for arbitrary /j, then gives the expression for ClHf\z,Q^/fi'^,Q^/M^). Presumably 
the ACOT scheme works at higher orders in a similar manner, with the higher order heavy quark 
coefficient functions being calculated explicitly (but needing explicit subtraction of divergences in 
(Q^/M^) beyond leading order). However, we note that ACOT do not usually use the scale choice 
/x^ = as we do. More common is fj? = + 0.5(^^(1 — M^/Q^)^ [23], which grows more slowly 
than our choice from the same value at = and is /x^ = O.SQ^ asymptotically. 

ACOT claim a smooth transition from the FFNS at low order by order. Their "LO" 
expression for the structure function is 

F^o^x, Q^) = Cl^'S'''{QyM') ® (i/(^2) + F(^2))o 

+ ''"%^;;^^'^ 4rgf^(QV/^^QVM^)0go^+^M^)■ 

^ (4.24) 

where from (3.11) and (4.23), 

Clnf'i^, QVi^',QVm') = Cl^f\z, Q^M') - {Hii'/M^) + c,,)P°, ® z5{xo -z){l + ^) . 

(4.25) 

There are a number of odd features associated with these expressions. Firstly, the "correct" 
threshold behaviour comes about only from a conspiracy of cancellation. Neither term in (4.24) 
respects the physical threshold individually and C2Hq''^{^^ Q'^ I ^ •, Q"^ /M"^) has a part with a thresh- 
old in and a part going like (4.23). In fact, since the first term in (4.24) grows more quickly 
than the subtraction term in the second term in (4.24), there will be nonzero (albeit very small) 
heavy quark structure function for W"^ < 4M^. Once all the necessary cancellation has taken place, 
the result is very good. This can be seen in fig. 8 of [23], and also in fig. 1 which is calculated using 
the ACOT "LO" prescription, our choice of renormalization scale, and the parton distributions 
obtained from our best fit (see later for details). There is a smoother transition in fig. 8 of [23] 
than in fig. 1 because their complicated choice of scale leads to /i^ departing slowly from and 
staying well below and hence to the growth of the charm parton distribution being effectively 
much slower than for the simple = jj? choice. The effect of the choice of renormalization scale 
on the speed of departure of the ACOT result from the LO FFNS result can be seen nicely in fig. 
1 of [26]. 

However, even though the cancellation of terms works well, (4.24) is at odds with the usual 
way of defining a LO expression, which usually only involves zeroth coefficient functions convoluted 
with the parton distributions obtained from the one-loop evolution equations. It is clearly of 
mixed order, and indeed, part of the expression is in fact renormalization scheme dependent, which 
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is certainly not correct for a LO expression. If we go to >> M^, (4.25) does not reduce to any 
fixed order expression in the ZM-VFNS. The first term in (4.25), represented by the dotted Hne in 
fig. 1, becomes the LO expression in the ZM-VFNS, but the second belongs to the NLO expression. 
One can see in fig. 1 that the total LO ACOT result is significantly different from the asymptotic 
ZM-VFNS result even at = lOOOGeV^. Similarly the derivative of (4.24) leads to terms both of 
C(Q!s,nj.+i) and C'(Q!s,„j,+i), and will have a renormalization scheme dependent part. This mixing 
of orders is not acceptable. 

Alternatively, with the choice of C^HH^''^i^^ Q^/M^) made, the usual way of ordering the ex- 
pansion for a structure function leads to serious problems. Using what one would normally consider 
the LO expression, F^jjix, Q^) = C^H^^'^iQ^ /M^)®{H{fi'^) +H{n'^)), has only a sharp threshold 
in and the rate of growth of F2^h{x, Q^) would be very discontinuous at = and a great 
deal too fast just above this. This can easily be seen on fig. 1 where this contribution is represented 
by the dotted line and labelled "charm quark" . It deviates very quickly from both the continuation 
of the FFNS expression and from the total expression. Using the NLO expression ordered in the 
usual manner^ the effect would be lessened, but would still be significant. The subtraction piece in 
^2Hg'^i^^ <5^//x^, Q^/M^) would largely cancel the quick growth generated by the LO evolution of 
the charm parton distribution, but the NLO evolution would still be uncancelled. This effect can 
be seen in fig. 9 in [23] where NLO parton distributions are combined with what is called the LO 
coefficient functions and in fig. 2, where we do the same thing using our parton distributions and 
jx"^ = Q"^. Here the subtraction term only partially cancels the charm quark contribution and the 
total quickly departs from the continuation of the FFNS structure function, and the effect increases 
at smaller x. The all-orders definition of the coefficient functions in (3.5)-(3.9), if indeed it is the 
all orders definition in the ACOT scheme, guarantees that the correct low Q"^ behaviour will be 
restored when working to all orders, but in this scheme this behaviour will come about only due 
to the mixing of effects at different orders. At low orders the discrepancy is still large. We note 
that the MRRS scheme [28], which incorporates mass effects into the evolution, but has a similar 
definition of coefficient functions to ACOT (though with the usual ordering), suffers badly from 
this problem outlined above. At the transition point where the heavy quark starts contributing to 
the heavy quark coefficient function directly there is a very distinct kink, and the total rises very 
quickly above the continuation of the FFNS expression, as seen in figs. 6 and 7 of their paper. 

We do not believe that the method used by ACOT (or MRRS) is a satisfactory way in which to 
define the coefficient functions in a VFNS, and we certainly do not believe that it is Tinique. It is a 
choice, as our prescription is a choice, and as we have discussed in §3, we do not believe that any are 
strictly "correct" . However, using the ACOT choice the calculation of the heavy quark coefficient 
functions proceeds as though the heavy quark parton distribution is due to intrinsic presence of 

^ It is an expression of this general form that is used in the recent global fits to data [27]. 
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the heavy quark rather than it being generated from (at least mainly) the gluon. In particular 
the heavy quark coefficient function contains no reference to the kinematic threshold in W^. This 
necessitates a mixing of orders to get satisfactory results. We believe it is far more useful to choose 
the heavy quark coefficient functions so that they reflect the physics and all automatically contain 
at least the correct form of low behaviour, and our prescription guarantees this. 



5. The VFNS in Practice. 



We now discuss how our procedure is implemented in practice. Of course, in practice the 
first heavy quark we encounter is the charm quark with rric ~ l.SGeV. First we consider the LO 
expression. Denoting e = ml/Q^ the LO FFNS heavy quark to gluon coefficient function is 



+ {8z{l - z) - 1 - Aezil - z))v 



(5.1) 



9{W^ - Ami) 



where W'^ = Q'^iXjz — 1), the gluon quark centre of mass energy, v is the velocity of the heavy 
quark or antiquark in the photon-gluon centre-of-mass frame, defined by = 1 — 4m^/Vr^, and 
Pqg{z) = {z^ + (1 — zY)- These u-dependent terms ensure that the coefficient function tends to zero 
smoothly as W'^ = 4m^ is approached from below, and hence the structure function has a smooth 
threshold in W"^. Taking the In Q^-derivative of this is a straightforward matter and results in 



dC2,^g^{z,e) 
dlnQ2 



(P,V-) + 2e^^-^-16eV)l 
1 — z V 



+ (-4ez(l - 3z) + 16e^z^) ln(^ - 
and it is easy to see that in the limit 



l + v 



+ (4e2;(l - z))v 



(5.2) 
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P° (z) 



(5.3) 



Hence, from (4.6), we see that C'^^'^'^'°(z, e) must indeed tend to the simple form z6{l — z) in this 
limit. 

Solving (4.6) for C^^'^'^''^(z, e) at arbitrary e is not too complicated. Taking moments of both 
sides the Mellin transformation of Co^J.^^'^{z, e) is the product of the Mellin transformation of 



^2,cc 



dC. 



^'j°^gii' and the inverse of the Mellin transformation of {P^g{z)), which is 

2 



/ 

^0 



z^-^P^^g{z)dz = 



n 



n 



n(n + l)(n + 2)' 



(5.4) 
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This leads to the following expression for the convolution of C'^^'^'^'°(^;, e) with the heavy quark 
distribution, 

Jx otlnc^ \ / dyx/z) 

where xq = (1 + 4e)~^ and r{z) is given by 



1 

r{z) = z2 



fV7. i\ 3 . /Vr, 1^ 

cos — — In - H — sm — — In - 

V 2 zJ ^ \ 2 z. 



(5.6) 



Using this expression we are able to calculate the LO contribution to the heavy quark structure 
function using a particular set of parton distributions. In practice we use those obtained from a 
global fit to structure function data using the NLO formalism (details later) . In order to get the LO 
parton distributions we simply take the same input parameterizations for the partons but evolve 
them using the LO evolution equations. Our prescription for the LO agiCf') across the charm 
threshold is to define 

and 

a,,4(Q') = a.(Q^4), (5.8) 
i.e. kqcD is defined for 4 flavours, and take for three flavours 

= «7'(Q',3) + a;'{mlA) - a;'{ml,3). (5.9) 

This prescription precisely reproduces the results of summing the leading logs in {Q"^ /m^) in (3.14). 
The results of the LO contribution for the heavy quark coefficient function are shown in fig. 3, along 
with the continuation of the LO FFNS expression and also the LO ZM-VFNS expression. One can 
see that the LO VFNS expression departs very smoothly from the continuation of the LO FFNS 
expression, then rises above it, and in the limit of very high becomes essentially identical to the 
LO ZM-VFNS result. This is precisely the behaviour we would expect. We also note that unlike 
other approaches the expression does not rely on any cancellation between terms. 

We now consider also the NLO expression for the heavy quark structure function. As well 
as the LO coefficient function just introduced we include the C(Q:s,n/+i(Q^)) coefficient functions. 
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The expression for C^^'^(z, e) is as in (4.7), and is in terms of quantities we have written exphcitly 
above in (5.1) and (5.2), i.e. 

C2,cg e) = C'2,cg ^) " HQ Mc) tiin(Q2) ' (^•^^) 

Hence, there are no new problems in implementing C^^'^(z, e). In the limit — ^ oo the first of 
these becomes 

Cl^,;\z,e)^P%{z)\n{^^^^^+{Sz{l-z)-\). (5.11) 
Using this, along with (5.3) and the definition (5.10), we see that in the limit e ^ oo 

Ci:,\z,e) ^ nV-)l-(^) + (8^(1 (5 ^2) 

in the MS renormalization scheme. 

From the arguments leading up to (4.12) it is clear that C^^'''^'^'^(z, e) also tends to the 
correct asymptotic limit, and indeed, all coefficient function are constructed so that this will be 
true. However, it is not possible to exhibit this in such an explicit manner since the expression 
for C;^^'^'^'^(z, e) depends on C^^'^(2;,e) for which there is no analytic expression.^ Likewise, it 
would be extremely difficult to implement C^cc^^'^{z, e) into the calculation precisely. In practice 
we find that the contribution to the total heavy quark structure function of this term convoluted 
with the heavy quark parton distribution is only a very small fraction of the total. Hence, we 
include this contribution to the total heavy quark structure function in an approximate manner, 
being confident that it is very far from being significant at the present level of accuracy required. 

Using our NLO prescription we use our NLO partons to calculate the heavy quark structure 
function. Our prescription for the NLO as{Q^) is to define as{Q^,n) by the standard two loop 
extension of (5.7), and then to use (5.8) and (5.9) once again. (5.9) does not sum all leading and 
next-to-leading logs in (Q^/m^) absolutely correctly, but is an extremely accurate representation 
of the precise expression. The NLO charm structure function is shown in fig. 4 along with the 
continuation of the NLO FFNS expression and the NLO ZM-VFNS result. As at LO the VFNS 
departs very smoothly from the continuation of the FFNS expression. Although at this order we 
have not been able to demand absolute continuity of the derivative of the structure function across 
= ml we see that there is no visible evidence of discontinuity at all. In fact the transition from 
one scheme to the other is smoother than at LO. Also the VFNS stays closer to the continuation 
of the FFNS at higher at this order. This is as we would expect, since as one works to higher 

^ We are grateful to Jack Smith and Steve Riemersma for providing the extensive program to compute 
the 0{a'i) FFNS coefficient functions [29]. 
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orders in in the FFNS one automatically includes more ln((5^/m^) terms which are completely 
summed in the VFNS. At all orders the two schemes become identical insofar as physical quantities 
are concerned. At very high then our expression tends towards the NLO ZM-VFNS exactly as 
required. 

Of course, at higher we also have effects due to the bottom quark which has m& ^ 5GeV. 
Below = ml there is no bottom quark distribution and we take account of the bottom quark 
effects by using treating all diagrams including bottom quarks in the fixed flavour scheme, and all 
other effects decouple. At = ml we switch to a VFNS for inclusion of the bottom quark effects. 
Analogously to the charm quark this involves switching to a 5 flavour coupling constant, defined 

by 

a-l{Q^)=a-\Q^5) + a-\ml4)-a-\ml5), (5.13) 

and beginning the evolution of the bottom quark density. The VFNS cocflicient functions are 
defined using a generalization of (3.5)-(3.9). There are now two heavy quarks in the definition of 

FF S FF S 

the fixed flavour number scheme, so we have two extra equations for C2 5^ and , and there 

\^ F S FF PS 

are now extra VFNS coefficients such as or C2 . Also, the finite operator matrix elements 

will depend on both the charm and bottom mass effects. However, exactly the same principles as 
outlined in the last section apply for determining the VFNS coeflficient functions. At low orders 
in as{Q^) there is no mixing between the charm mass effects and the bottom mass effects. Hence, 
the VFNS charm coefficients functions we have mentioned explicitly above remain the same above 
= ml (except for a completely negligible change in C^^'^*^'^) and the bottom quark coefficient 
functions are obtained simply by replacing mc with m^ and n/ + 1-ff avour splitting functions with 
n/ + 2-ff avour splitting functions. At higher orders the VFNS charm coefficient functions change 
above Q"^ = ml, acquiring m^ dependent corrections (in particular C^^^^^ comes into existence), 
and bottom coefficient functions acquire charm mass corrections. 

Thus, we have described how one may implement our prescription for the VFNS in practice, 
showing that there is no real difficulty. We have also demonstrated that the results have precisely 
the properties that our theoretical arguments in the previous section lead us to expect. In order to 
make even more concrete statements regarding the suitability of our VFNS for the calculation of 
structure functions we will now discuss a comparison with data. 

6. Phenomenological Results. 

Using the prescription for heavy partons discussed above we can calculate the full structure 
functions in terms of input parton densities for the light quarks and gluon. The input scale is chosen 
as Qo = IGeV^, and the input parton distributions are then determined by performing a best fit 
to a wide variety of structure function and related data. Hence we repeat the type of procedure 
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adopted by MRS and CTEQ (and others) in their global determination of parton distributions. We 
note that apart from the masses rric and nib, which we fix at 1.35GeV and 4.3GeV respectively, the 
heavy quark contributions to structure functions have no free parameters. The overall description 
of the data is shown in table 27 There is a clear improvement when compared to the MRS(R2) 
fit [30], which uses the ZM-VFNS prescription with a phenomenologically motivated smoothing 
function, and a small but definite improvement when compared to the MRRS fit.^ 

Pull details of a global analysis using this procedure will be presented in a future publication. 
However, here we will concentrate on those data which describe the charm contribution to the 
structure function only. The component F2,c has been measured at intermediate x values by 
EMC [32] (via the detection of inclusive muons) and at small x by the HI [4] and ZEUS [5] 
collaborations at HERA (through measuring D and D* cross -sections). The latter indicate that 
as much as 20 — 25% of the total F2{x, Q'^) is due to F2^c{x, Q^)- While these data on F2,c(a;, Q^) 
are not included as part of the general fitting procedure we can compare them with our theoretical 
predictions. The results are shown in fig. 5. A very good description of both the small and medium 
X data is achieved for a charm mass of rric = 1.35GeV,^ although there is a strong suggestion that 
a slightly higher mass would be preferred, i.e. the curves for rric = l-SGeV give a rather better 
description. In fact it is the data for ~ 2GeV^ which strongly favour this higher value of rric. 
Since in this region of the theoretical approach is unambiguous, i.e. the true result must be 
very similar to the FFNS calculation, these points may be thought of as determining tlic value of 
rric at about 1.5GeV, and a value very similar to this that should be used over the whole range. 
Also, in fig. 6 we show the relative importance of the charm and bottom components to the total 
structure function, and note the large fraction which is due to charm in the HERA low x region. 
The bottom contributes no more than ^ 4% in any currently accessible range of x and Q^. 

7. The Longitudinal Structure Function. 

Finally we discuss our prediction for the heavy quark contribution to the longitudinal structure 
function. Although there is no data directly available on this quantity, we feel that it is an important 
issue. This is because the implementation is a little different from the case oiF2{x, Q^), also because 

We note that we do not alter the values of F2(x,Q^) for the HERA data to take account of our 
predictions for Fl{x,Q^), as should really be done. The Fl{x,Q^) values used in [2] and [3] are obtained 
using a NLO-in-as((5^) calculation, and so are not very different from ours in general, and the number of 
points affected is relatively small. Hence the quality of the overall fit is very insensitive to the neglect of 
this small correction. 

* The fit is not as good as the LORSC fit [31] [20], which includes ln(l/x) corrections but not yet charm 
mass corrections. 

^ There is also a single EMC data point at x = 0.422 and = 78GeV^ not shown in fig. 5, which has 
F2,cix, Q^) = 0.00274 ± 0.00152 compared to a prediction of 0.0003. 
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the charm contribution has a very important bearing on the total longitudinal structure function, 
contributing up to about 35% in the ZM-VFNS for example, and finally because our results are 
very different from those in other ways of implementing a variable flavour number scheme. 

As for ^2(0;, Q^), and for the same reasons, there is no way we can obtain the (hypothetical) 
absolutely correctly ordered expression. Therefore, as in this previous case we want some relatively 
simple prescription which will reflect the physics involved correctly. There is a lot of similarity 
between our order-by-order prescription for the heavy quark contribution to Fl{x,Q^) and that 
for F2{x,Q^), and the equations that the VFNS coefficient functions must satisfy in terms of the 
operator matrix elements and the FFNS coefficient functions are once again (3.5)-(3.9). One large 
difference between the two is the fact that in a zero-mass formulation there are no zeroth order 
in as{Q^) coefficient functions for the longitudinal structure function, and hence the 0{as{Q^)) 
coefficient functions are leading order and renormalization scheme independent. All previous imple- 
mentations of a VFNS [22] [23] [28] have included a zeroth order heavy quark longitudinal coefficient 
function, i.e. the term in (4.23) which oc 4M^/(5^. This procedure means that there is a coefficient 
function at lower order than the one which becomes leading order in the ZM-VFNS limit, and 
hence in order to reach this limit with the LO VFNS expression coefficient functions at both zeroth 
and first order in as^uf+iiQ'^) would need to be included. Also, if one includes any zeroth order 
coefficient function, using the expression (3.5) for the O^as^nf+iiQ"^)) gluon coefficient function 
results in C^^]^{z, /M^) having a component which is renormalization scheme dependent. 

Hence, wc choose not to have any zeroth order contribution to the longitudinal coefficient 
functions. As with F2{x,Q'^) our VFNS coefficient functions are then determined entirely by the 
requirements of reduction to the ZM-VFNS order by order as ^ 00 and continuity with the 
FFNS across the boundary = M^. Therefore, the prescription for the VFNS longitudinal 
structure function is very similar to that for F2{x,Q^), except that the relative order of heavy 
quark coefficient functions above = and light quark coefficients at all is one higher, i.e. 
0{as,nf+iiQ^)) is leading order etc.. The prescription for the FFNS structure function at fixed 
order is then very straightforward, i.e. 

Fs,M') = EEi^''^) c^r"'^\M-m®f:ii{Q% n=i^^, (7.1) 

m=0 a ^ ^ 

for both the heavy and light quark structure functions. The general form of the expression above 
= is the same as this, i.e. 

FU-^Q') = ) Clfr-'+'iMyq^) /^^+\Q^), n = 1 ^ 00. 

m=0 a ^ ^ 

(7.2) 

Since the expressions are now of an identical form both sides of the transition point (which was 
impossible for F2{x, Q^) because of the requirement of zeroth order heavy quark coefficient functions 
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above Q"^ = M^), and are identical to all orders, continuity of the structure functions themselves 
is guaranteed order by order in as{Q^) across the transition point. However, as in the case of 
F2{x,Q'^) the heavy quark coefficient functions at each order have to be determined using some 
prescription. (This ambiguity has no effect on the continuity of the structure function since at nth 
order in the expression for the structure function the n^^j-order heavy quark coefficient functions 
only appear coupling to the zeroth order heavy quark distribution, which vanishes at = M^.) 

As with F2{x,Q^) it would be nice to demand both continuity of the structure function and 
its In Q^-derivative across the transition point. Since the expressions for the structure function are 
of exactly the same form both above and below = in this case (essentially because there are 
no zeroth order terms in the longitudinal structure function), we can now attempt to equate the 
In derivatives of the nth order terms in both expressions rather than making the derivative of 
the ntfe-order VFNS expression match on to the nth order derivative in the FFNS as was necessary 
for F2{x,Q'^). As in the previous case wc have enough choice to demand only that this is true in 
the gluon sector, but again by far the dominant contribution to this derivative comes from this 
sector. However, now we have an additional problem compared to the previous case. This can be 
seen by examining the lowest order expressions. 

In the FFNS the lowest order expression for the heavy quark structure function is 



97 iQ' 



(7.3) 



while in the VFNS it is 



2x _ (^s,nf+l{Q'^) ^VF,1 1 



VF,1 



2tt 



and from (3.5) we have 



(7.4) 
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(7.5) 



for the charm quark, where again e = m^/Q'^ and 9(W'^ — 4m^) is implied whenever v appears. 
Clearly the structure functions arc the same at = M^. It is also clear that the 0(as{Q^)) 
expression for the In Q^-derivative is the same on both sides, i.e. 



dlnQ2 



27r 



dlnQ2 



(7.6) 



However, this expression, which vanishes as Q^/M'^ — > oo, is lower order than the leading order 
asymptotic expansion, which is 0{al{Q^)). It is this, rather than a zeroth order coefficient function, 
which truly reflects the fact that the heavy quark longitudinal structure function has behaviour 
which begins at lower order than the massless expression. This 0{as{Q^)) derivative means that 
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while the asymptotic 0{a1{Q'^)) expression for the total derivative of the longitudinal structure 
function is renormalization scheme independent, it contains a part which vanishes as Q^/M^ — oo 
which is renormalization scheme dependent. This is different to the case for F2^h{x,Q^) where 
the leading asymptotic and 0{M^/Q'^) parts of the derivative are both of the same order, i.e. 

If we treat the 0{as{Q'^)) component of '^^'^^"^q? ^ as a super leading part which is trivially 
continuous across = M^, and then examine the form of the 0(ag(Q^)) terms coming from the 
derivatives of (7.3) and (7.4), then since each of the leading-order expressions is renormalization 
scheme independent then so are the contributions obtained in the expressions for the derivatives. 
ExpHcitly we obtain in the FFNS 



and in the VFNS 



dlnQ2 \^ 27r 

+ P^,^K'iQ'))+crHHiQ'/M')0{p^,0g7^\Q') 

+ ® (iJ(Q2) + HiQ^)),) - l37+\HiQ') + H{Q'))o . 

(7.8) 

Prom previous arguments it is clear that the terms oc C^^^g{z, Q^/M^) in each equation are equal at 
= M^. The vanishing of the heavy quark distribution at this scale, lead to the single condition 

CrS(^>l) = 0, (7.9) 

in order to match these C(a^(Q^)) contributions to the derivative. Thus, we have this condition, 
along with the fact that C^^]j{z, Q'^ /M'^) must reduce to the correct asymptotic form, in order 
to determine C^^l^iz, Q^/M^). It is clearly possible to choose forms for C^^]i{z, Q^/M"^) which 
satisfy these conditions, but there is rather less guidance as to the precise form required than for 
P2,h{x, <5^), where the condition at = contained a component which was clearly identifiable 
as the asymptotic expression. 

This indeterminacy is due to the fact that the 0{M'^/Q'^) contributions to the derivative 
begin at one lower order than the asymptotic form, rather than our chosen manner of imposing 
the matching. If we had chosen to match (7.8) to the total 0{a1{Q'^)) expression for "^^^^"^q? ^ 
in the FFNS, rather than just the part coming from (7.7), i.e. analogously to F2^h{x,Q'^), we 
would have encountered a different problem. In this case the determined value of C^^\j{z,l) 
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would have contained a part oc — ^'"din'g'^ ' 'which contains C^^'^'^~^^'^ (z), and the asymptotic 

limit would therefore appear more naturally. However, the full expression for — ^'"dinQ^ ' ^^'^ 

consequently the full expression for C^^];^(z, Q^/M^) implied, contains 0{M^ /Q^) parts which 
are renormalization scheme dependent (since they are subleading to the ©(^^(Q^)) expression). 
This is not satisfactory in the definition of the leading order VFNS coefficient function, and the 
renormalization scheme dependent part of the expression should be removed. However, there is no 
unique way to do this, and hence the definition of C^^\j{z, M"^ /Q^) would be just as ambiguous 
as when using our chosen matching condition. 

Hence, we have to live with the fact that there is no completely satisfactory way to determine 
C]^^^(^, Q^/M^) from physical arguments. We choose to impose (7.9), as well as the fact that 
C^^\j{z,Q'^ /AP) must reduce to the correct asymptotic form, and also choose the coefficient 
function so that a smooth threshold in W"^ is automatically incorporated. A simple choice satisfying 
all these requirements is 

CrHHiz, QVM^) = ^v{l - MyQ^)z. (7.10) 

In practice this ambiguity has little effect phenomcnologically since the vast majority of the LO 
expression for ^{x,Q^) comes from the gluon contribution which is determined uniquely. Using 
(3.5) we have now also defined C^^^(2;, Q^/M^), i.e. 

<5^(., QVm^) = Cl^^liz, Q'/M') - HQ'/M')P^, <^^(QVM^), (7.11) 

in MS scheme, although we do not have to make use of this yet. The fact that C^^\j{z,Q^ /M^) 
reduces to the correct asymptotic limit guarantees that C]^^^(z, Q^/M^) does. 

As far as the light quark contribution is concerned the coefficient functions are identical above 
and below = M^, and the only effect is the change of the evolution of the parton distributions 
and in the running of the coupling. The In Q^-derivatives of these LO light quark distributions 
which are entirely of 0{a^{Q^)), are not quite continuous across the transition point because of 
the flavour dependence of (3q and of the lowest order splitting functions, i.e. of i-*^g(2;). As in the 
0{a'^{Q^)) derivative for F2^i{x,Q'^) there is continuity in the gluon sector, but not in the quark 
sector. Phenomcnologically the discontinuity is very small, and becomes formally smaller as we 
work to higher orders. 

The result of our leading order calculation of F^^cix, Q^), using the same LO parton distribu- 
tions as before, is shown in fig. 7, along with the LO FFNS and the LO ZM-VFNS results. As in 
the case of F2^c{x, Q^) one can see that the transition from the FFNS result is extremely smooth, 
and of course, the the correct asymptotic limit is reached. We note that at low the VFNS result 
for Fl^c{x, Q^) is very different indeed from that in the ZM-VFNS. This leads to a very significant 
difference between the results for the total Fl{x,Q^) in the two different schemes, and important 
phenomenological implications. We also show explicitly the contribution made by the charm quark 
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distribution. At high this is unambiguously defined, and at low it is very small indeed. 
Therefore, the ambiguity in the low heavy quark contribution is not really significant. 

A very important difference should be noted between this approach and previous VFNS ap- 
proaches at this point. As already mentioned, all previous approaches have used a zeroth order 
charm quark coefficient function of the form 

C'K°(^,6)=46Z(5(io-^). (7.12) 

If one were to regard the LO expression for Fl^cIx, Q^) as just this coefficient function convoluted 
with the heavy quark distribution function then the behaviour would be rather strange, having a 
sharp threshold at = ml, growing quickly, then turning over and going to zero as Q^/m^ — > oo. 
If, as is more likely, the LO expression is taken to include both the zeroth order and 0(as{Q'^)) 
coefficient functions, so that the correct asymptotic LO limit is reached, then C]^^g{z, e) is defined 
by (3.11), i.e. 

= C^,\z, e) - ® z5{xo - z){\n{Q^/ml) + c.,). (7.13) 

As well as this introducing incorrect renormalization scheme dependence into a leading order ex- 
pression (via c^s), it has unfortunate phenomenological consequences. The VFNS differs from the 
FFNS expression by 

Aez6{xo -z)® [(c(Q2) + c{Q^))^ - a.,4(<3')n°9 HQ^/ml) ® gt{Q^)] , (7.14) 

where we have used MS scheme. These two terms are intended to largely cancel at and just 
above Q"^ = m^, ensuring a relatively smooth transition as in the ACOT prescription for the LO 
expression for F2^c{x,Q^)- The procedure works well in the case of F2^c{x,Q^), and the transition 
is quite smooth, as we have seen. However, the cancellation is not exact (otherwise we would just 
have the FFNS), (c(z,Q2) + ciz,Q^))o « as,iiQ^)P!lgHQyml) g^iQ^) for just above ml, 
but the resummation of the logs in the evolution of the charm quark distribution leads to differences 
appearing. To a rough approximation. 



{c{z, Q^) + c{z, Q^))o - a,AQ^)P° HQVrnl) ® g^{Q^) 



(7.15) 



at moderate Q"^ . Inserting this into (7.14) leads to 



4e(a,,4(g') HQ^/ml)fz5{xo - z) ® P^ ® ® gt{Q^). (7.16) 

For 5 — lOGeV^ this expression is comparable in size to the FFNS component of the full 

expression for FlAx,Q'^), which ~ asg^iQ^) with damping due to kinematic factors (and which 
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is more than 10 times smaller than the LO FFNS component for F2^c{x,Q^))- However, it falls 
away quickly at larger Q^. This leads to the LO VFNS expression for Fl c{x,Q^) increasing very 
quickly above the FFNS expression above the transition point, dramatically slowing, or perhaps 
even falling at ^ 5m^, and then smoothly approaching the correct asymptotic limit, i.e. there 
is a very pronounced unphysical bulge in the value of Fl^c{x, Q^) calculated in this way. When one 
calculates Rc = i^L,c/(-^2,c — -^l,c), which exhibits the relative rate of growth of Fl^c and ^2,0 the 
effect is demonstrated much more clearly as a distinct hump peaking at about = 3m^. This can 
be seen very clearly in fig. 9 of [28], where the effect is particularly dramatic since the evolution 
of the heavy quark distribution there is even quicker than in MS, and is at NLO. However, the 
treatment of coefficient functions follows the same general principles as ACOT, and the same type 
of effect, if somewhat smaller (the reduction depending very strongly on the particular choice of 
renormalization scale - one which departs extremely slowly from as increases could remove 
the effect) will be clearly seen in their expressions.^^ Even in the absence of detailed data this type 
of effect seems sufficient to rule out this approach as a suitable way to order a VFNS expression. 

We now consider the NLO expressions for the longitudinal structure functions. For both 
heavy and light quark structure functions both above and below = we add to the LO 
expressions the 0{a'l{Q^)) coefficient functions convoluted with the LO parton distributions and 
the 0{as{Q^)) coefficient functions convoluted with the NLO parton distributions. Let us first 
consider the heavy quark coefficient function. It is guaranteed by satisfying (3.5)-(3.9) order by 
order in as(Q^), while also satisfying the correct relations between parton distributions and the 
coupling, that this procedure will lead to structure functions which are continuous across = . 
This is straightforward, if a little tedious to check. Continuity of the derivative of the heavy quark 
structure function across the threshold is not guaranteed, but depends on the particular choice of 
the heavy quark coefficient functions. We can compare the derivatives of the full NLO expressions 
in both the FFNS and VFNS up to 0{al{Q^)). From the conditions we have already imposed using 
(3.5)-(3.9) it is guaranteed that all new terms we introduce which behave like a^{Q^), i.e. those 

,^F(V)F,2 

depending on rfinQ^ > ^^^^ continuous across the transition point. Again this is straightforward 
to check. If we examine the 0{a^{Q^)) contributions to the expressions then both in the FFNS and 
VFNS these are very involved, i.e. containing rather more terms than (4.10) and (4.11). However, 
as with (4.10) and (4.11) many of these terms vanish at = M^, because the heavy parton 
distribution vanishes here, also because in this case C]^^^^'^{z,l) = 0, and also because many 
other terms are the same in both expressions. A long, but entirely straightforward calculation 

In fact, since as we see in fig. 1 at a; = 0.005 the subtraction term is larger than the heavy parton 
distribution, the effect will be negative. 
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reveals that if we require continuity of the derivative in the gluon sector we have the requirement 
that 

(7.17) 

at = M^. So at 0(a^(Q^)), as at 0{as{Q^)), there is no impUcation of the asymptotic form 
required of the heavy quark coefficient function in the condition at = , but the condition is 
no longer that the coefficient function is zero at this value of Q"^. We can understand where the 
nonzero terms come from quite easily. If we had used the whole of the 0{al^n^{Q'^)) expression 
for the derivative of the heavy quark structure function in the FFNS, and equated this to our 
VFNS expression then the asymptotic form of C^^^jj{z, Q'^/M'^) would have appeared naturally in 

the expression for ^^^q! ■ However, by examination of the expression for ^^^q'! contained within 
(3.5) we would find that the definition of C]^^^^''^{z,Q'^ /M"^) would also need to contain terms 
of the sort in (7.17), as well as others which vanish at = M^, in order to reduce to the correct 
asymptotic limit. However, in an analogous fashion to our previous discussion at leading order, we 
do not use this technique since parts of the 0(M^/Q^) corrections to ^^^gi are properly of NNLO, 
i.e. are renormalization scheme dependent in such a way as to compensate for the renormalization 
scheme variation of the NLO terms. This would require an ambiguous subtraction procedure for 
these terms, and we would have no more real information than that contained in (7.17) and the 
asymptotic condition. 

Hence, as for the 0{as,nf+i{Q'^)) coefficient function we make a simple choice for the coefficient 
function which satisfies (7.17), which reduces to the correct asymptotic limit, and which explicitly 
contains the correct threshold behaviour. Once again we multiply the asymptotic limit, which makes 
no appearance at = M^, by {1—M^/Q'^)v. We multiply the terms appearing in (7.17), but which 
must disappear asymptotically, by M'^/Q'^ (in this case the threshold behaviour is automatically 
contained in the expressions). Hence we obtain 

crHrMVM-) =(i - ^)vc2^,:;'\z) + ^{p^,r ® [ci'hUqvm-m^^' - ^-^) 

(7.18) 

This definition is ambiguous at low Q^, but as at leading order the total heavy quark structure 
function at NLO is totally dominated by the gluon contribution. We also note that the ambiguity 
introduced at LO from the definition of the heavy quark coefficient function is very largely negated 
at NLO by the inclusion of Cl^^^^'^{z,Q'^ /W^) in the expression for C^^'2(z, Q^/M^) (7.11). As 
we work to higher orders the ambiguity formally disappears. We also note that the coefficient 
function C]^5&^ 

(z, Q^/M^) is the sum of the nonsinglet and pure singlet coefficient functions. We 
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are free to separate them as we wish, using the condition that each tends to the correct asymptotic 
hmit. It would also be desirable to choose each so that they respect the kinematic threshold. The 
choice has no bearing on the expression for the structure function but a simple choice is to let the 
nonsinglet part contain all parts oc v, and to split the other part simply in terms of the asymptotic 
form. 

Comparing the 0{a^{Q'^)) expressions for the derivative of the heavy quark structure func- 
tion which are proportional to the singlet quark distribution, then as for the NLO derivative for 
P2,h{x,Q^) we see that continuity is not achieved. The difference between the VFNS expression 
and the FFNS expression is 

2(/3o^ - Po'^')Cl57'"(Q"/M") ® ^o^(Q')' (7-19) 

where C^f ^^'^(z, QV^^) = C'lSf ^'^(^' <5V-^^)- This NLO effect is very small, and as for 
P2,h{x,Q^) the effect disappears as we work to higher orders. 

For the light quark structure functions there is one choice to make. There is a mass de- 
pendent contribution to the nonsinglet coefficient function at 0{al{Q'^)), but the form of the 
VFNS coefficient function is determined entirely by (3.7). In essence the mass dependent cor- 
rection to C^^''^f''^{z) contain a piece which becomes constant asymptotically which represents 
(J NS,nf +1,2^^-^ — C^^''^f''^{z), and a piece which grows like ln((5^/M^) which takes account of the 
difference between the n/ + 1 and n/ flavour couplings. The 0{al{Q^)) InQ'^ derivative is slightly 
discontinuous at = M^, but this is corrected by inclusion of the 0{al{Q'^)) coefficient functions. 
For the pure singlet and gluon coefficient functions coupling to light partons there are no mass de- 
pendent corrections in the FFNS at 0{a^ nj{Q'^))^ '^i^d wc simply use the same coefficient functions 
above and below = M^. Continuity of the NLO structure functions is then automatic. How- 
ever, the pure singlet coefficient function C^^]^'^'^(z, Q^/M^) becomes nonzero at this order. It 
can be determined by demanding continuity of the In Q^-derivative of the structure function in the 
gluon sector. This results in a similar procedure as for the heavy to heavy coefficient function: the 
asymptotic form is put in by hand and multiplied by (1 — M"^ /Q'^)v while there are nonzero terms 
at Q2 = which arc multiplied by v to ensure that they vanish as Q"^ /M"^ oo. The calculation 
is straightforward, and we do not present details here. As for the heavy quark structure function 
the C'(ag((5^)) InQ^ derivative in the singlet quark sector is slightly discontinuous at = Af^, 
but again this is corrected at next order by inclusion of the 0{q!1{Q'^)) coefficient functions. 

Now that the NLO prescription for the longitudinal structure function is completely defined we 
can examine the results. Using our NLO coefficient functions (7.18) and (7.11) and the NLO partons 
obtained from the best fit we calculate the NLO charm quaxk longitudinal structure function. This 
is shown in fig. 8 along with the continuation of the NLO FFNS expression and the NLO ZM- 
VFNS result. Once again the VFNS increases above the FFNS result very smoothly despite the 
discontinuity in the In derivative in the singlet quark sector which is now demonstrably minute. 
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At very high our expression tends towards the ZM-VFNS as required, but as at LO the two 
become very different at lower Q^. As in the case of F2(x, Q^) at NLO the difference between the 
VFNS and the continuation of the FFNS is reduced compared to the difference at LO for the same 
reasons. Hence we have every reason to consider our prescription for the longitudinal structure 
functions quite satisfactory. 

In fact we can compare to some data. At x < 0.1 the VFNS, ZM-VFNS and FFNS values 
for the total Fl{x,Q^) are very similar. However, the NMC collaboration have produced data for 
0.11 > X > 0.0045 and l.SGeV^ < < 20.6GeV^ [33], increasing as x increases. These data 
are in the region where our VFNS prescription produces very different results to the ZM-VFNS 
(but almost identical to the FFNS) for Fl^c{x,Q'^), and hence significantly different results for the 
total longitudinal structure function. Using the parton distributions obtained from our best global 
fit we produce predictions for R{x,Q^) = Fl{x,Q'^)/{F2{x,Q'^) - Fl{x,Q^)) using the ZM-VFNS 
and the VFNS, and compare to the data. The results are shown in fig. 9. The kink in both curves 
at the lowest x values comes about because for all data points other than that at x = 0.0045, as 
X decreases also decreases smoothly, while for this point the extraction of R{x, Q^) relies on an 
extrapolation and the is actually almost identical to that for the x = 0.008 point. The kink in 
the ZM-VFNS expression at Q"^ = = l.SGeV^ is due to all charm coefficient functions turning 
on discontinuously at this point. Comparing to the data it is clear that the ZM-VFNS curve is 
much too large for most of the x-range, while down to x = 0.01 the VFNS curve matches the 
data quite well. Thus, there is strong evidence for including charm mass effects in the longitudinal 
structure function, and our prescription seems reasonably successful. Other VFNS prescriptions 
would lead to R{x, Q"^) somewhere between the two curves. The theory is clearly below the data for 
the lowest x and points where the charm contribution to i?(x, Q^) is extremely small, i.e. the 
VFNS and ZM-VFNS curves are almost identical. The smallness of i?(x, Q^), and the decrease with 
decreasing x at constant in this region are largely due to a negative small x contribution from 
C2^-g{z) which becomes increasingly important as x and fall. Thus, the difference between the 
theory and data for the two lowest x points is perhaps a sign of the failure of the NLO-in-Q;s(Q^) 
calculation of structure functions at small x.^^ 

As with F2{x, Q^) the NLO calculation is the best that can be done explicitly with the present 
knowledge of structure functions. However, as in this previous case, we outline the procedure for all 
orders. The general form of the expressions is presented in (7.1) and (7.2), and for the heavy quark 

The curve labelled Rqcd in fig. 10 of the NMC paper [33] contains little information. For Fl{x, Q^) 
it uses a LO formula [34] (and hence does not contain the important NLO small x effect) which assumes 4 
massless quarks at all , along with a gluon which has been extracted using a NLO fit in the FFNS [2]. 
Moreover, this gluon is not constrained at large x and is highly inconsistent with large x data. Prom the 
momentum sum rule this means its form at small x is also much different to a well constrained gluon. 



36 



structure function nothing essentially new compared to the LO and NLO prescriptions occurs. At 
rifh nontrivial order we determine C^^^iz, Q"^ /M"^) by demanding continuity of the In derivative 
of the structure function in the gluon sector, and by demanding the correct asymptotic form. At 
each order the correct asymptotic form will not appear in the continuity conditions and need to be 
introduced by hand. Each time we multiply by {1 — M'^/Q'^)v. At each order there will also be terms 
introduced by the continuity demand which must vanish as oo, and we multiply these by 

M'^/Q'^. At every order this determination of C^^^^^iz, Q^/M^) predetermines C^^'^^'^iz, Q^/M'^) 
and C^^^'^^ {z, Q'^ /M'^) by using (3.5) and (3.6) to 0{a'^~^^{Q'^)). The comments concerning the 
separation of C^^^{z,Q'^ /M'^) into nonsinglct and pure singlet parts in §4 apply again. 

For the light quark structure function the procedure at higher orders is also straightforward. 
As with the 0{a'^,^^_^_-i^{Q'^)) coefficient function we determine C^^'^{z,Q'^ /M"^) by demanding 
continuity of the derivative of the light structure function in the gluon sector, analogously to 
the heavy quark sector. With this one degree of freedom eliminated in this way all other VFNS 
coefficient functions are determined uniquely order by order in ag^uf+i 

(Q2) by (3.7)-(3.9), i.e. this 

determination of Cl^';^{z,Q^ /M^) predetermines Cl^g''+\z,Q^ /M^) and C]^^^^'"+'(z, QV-^^) 
by using (3.8) and (3.9) to 0«+1(Q2)). 

Thus, we have completely defined our prescription for calculating the structure function 
Fl{x,Q^) order by order. As for F2{x,Q^) we can sum it up in a simple diagram, shown in 
table 3. The generalization to the case of two heavy quarks follows the same lines as for the case 
of F2{x,Q^) which was discussed at the end of §5. For < ml the bottom quark effects are 
all treated via FFNS coefficient functions while in the region above = ml wc have a variable 
flavour scheme for both the charm and bottom quark. For high orders in 

(Q^) there will be 

mixing of the effects of the two quarks, but for the orders currently available in practice the mixing 
is extremely small indeed, as with F2{x,Q'^), and the bottom coefficient functions are essentially 
the same as those for charm with m,c mi, and with 5 flavours rather than 4. 

Our prescription uniquely determines all VFNS coefficient functions, and as for F2(x,Q'^), 
while not leading to absolute correctly ordered expressions it is a relatively simple prescription 
for obtaining order by order structure functions which are very similar to the hypothetical strictly 
correct ones, which reduce to the correct asymptotic form order by order in as^uf+iiQ"^), and which 
are consistent with physical requirements order by order. All ways of satisfying both (3.5)-(3.9) 
and the correct asymptotic limits will be correct in a certain sense (provided they are consistent 
with ordering within a given renormalization scheme), but many will have behaviour which is 
unsatisfactory for not much larger than M^, and we have seen an example of this. As with 
F2{x,Q^) we believe our prescription to be very suitable. 
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8. Summary and Conclusion. 

In this paper we have constructed an order by order in Ug prescription for calculating the 
neutral current structure function including the effects of a massive quark. For the region < 
this has essentially just been the normal FFNS, where the heavy quark is not treated as a constituent 
of the hadron, but all heavy quarks in the final state are generated via the electroweak boson 
interacting with light partons. For > we have to solve the problem of summing large logs 
in Q^/M"^ and jU^/M^ which appear at all orders in as{n^). The easiest way to do this is to to 
treat the heavy quark as a parton, in which case the logs will be summed automatically when 
one solves the evolution equations for the partons. If one chooses the parton distributions above 
n"^ = to evolve as though massless and in the MS scheme, then the new n/ + 1 flavour parton 
distributions are determined in terms of the FFNS parton distributions at all by well-defined, 
calculable matrix elements which contain logs in fi^/M^. In particular the heavy quark distribution 
is determined entirely in terms of the light parton distributions. The matrix elements can then be 
used to define the n/ + 1 flavour parton distributions in terms of the n/ flavour distributions at 
some scale (in practice = is by far the most convenient), and the evolution upwards can 
take place in terms of n/ + 1 massless flavours with the correct asymptotic limits being guaranteed. 
If the massless n/ + 1 flavour coefficients functions are used then the correct asymptotic limit for 
the structure functions is also reached. 

The main problem lies in obtaining the correct description in the region not too far above 
= . We have demonstrated that this is achieved to all orders by defining the mass-dependent 
coefficient functions above = in terms of the operator matrix elements and the FFNS 
coefficient functions as in (3.5)-(3.9). However, we have also demonstrated that since there are 
more degrees of freedom on the right hand side of these equations than on the left, the additional 
ones all being coefficient functions coupling to heavy quarks, there is freedom in precisely how the 
coefficient functions may be chosen. Although in a true well-ordered calculation this ambiguity 
disappears, this manner of ordering is at the very least extremely complicated, involving parts of 
the FFNS from all orders in as(;U^) at each order in the calculation, and in practice is probably 
impossible, there being no unique prescription for ordering the ©(M^/Q^) terms. Hence we choose 
to order our calculation as in the normal ordcr-by-ordcr in ^^(/i^) manner, choosing the very simple 
natural scale fi'^ = Q'^, which puts all of the mass effects into the coefficient functions and guarantees 
the correct asymptotic limit order by order in UsiQ^)- We then determine the precise form of our 
heavy quark coefficient functions by demanding continuity not only of the structure functions 
at Q'^ = (which is automatic), but also the continuity of the In Q^-derivative of the structure 
function. In practice this exact continuity is only possible for those terms proportional to the gluon, 
but this is by far the dominant contribution. Our constraint then determines our prescription for 
dealing with heavy quarks completely, and incorporates the correct qualitative threshold behaviour 
into every coefficient function at each order of as{Q^), not relying on cancellations between terms 
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with incorrect behavior and of different orders to obtain satisfactory results. In practice the most 
important of our results are the zeroth order coefficient function for F2^c{x, Q^), (4.6), which exhibits 
the correct threshold behaviour in as well as reducing to the correct asymptotic form, and the 
absence of a zeroth order coefficient function for Fl,c{x,Q^), the 0{as{Q^)) coefficient functions 
being (7.5) and (7.10), which again exhibit the correct threshold behaviour and asymptotic limits. 

We display the results obtained using our prescription for neutral current structure functions 
in figs. 3, 4, 7 and 8, finding that they exhibit exactly the type of behaviour we would expect, 
i.e. smoothly deviating from the FFNS at low Q^, and tending towards the n/ + 1 massless results 
at high Q^, in all cases. In particular we notice that the bump in the charm quark longitudinal 
structure function at ^ lOGeV^ which occurs in other variable flavour number schemes is 
absent here. We also see that our predictions agree very well with the current data on the charm 
structure function which exists from 1.5GeV^ < < lOOGeV^, implying a charm quark mass of 
~ 1.45GeV. We note that comparisons of theoretical predictions with the complete range of data 
on the charm structure function appear very rarely (in particular, detailed comparison with EMC 
data is frequently omitted), and we strongly encourage this as the best constraint on any theory. 

The general technique can be applied to all other quantities in perturbative QCD which require 
the convolution of coefficient functions with parton distributions. We can always choose the parton 
distributions to evolve as though there are nj + 1 massless flavours in the MS scheme, factor 
these into the mass dependent operator matrix elements and the FFNS parton distributions and 
then obtain the coefficient functions in the variable flavour scheme in terms of those in the fixed 
flavour scheme by equating the parts proportional to each FFNS parton distribution. Indeed, 
the expressions (3.5)-(3.9) are not exclusive to neutral current structure functions, but apply to 
all quantities which can be written as the sum of convolutions of coefficient functions with single 
parton distributions. In the appendix wc discuss the case of the charged current structure functions 
as an example. For expressions involving more than one parton distribution the generalization is 
clear, e.g. for proton-proton scattering the FFNS and VFNS coefficients are related by equations 
of the form = Cy.^ A'^"' A'^^ . In all cases there will be ambiguity in definitions of the heavy 

parton coefficient functions, but these can always be eliminated by demanding as much continuity 
of the In Q^-derivative order-by-order in asiQ"^) as possible. 

Let us briefly discuss problems which arise in other approaches to heavy quark structure 
functions. Buza et al. do not provide a detailed prescription for the region of just above M^. 
They have a means of extrapolating the structure function from the FFNS result at < to 
the ZM-VFNS result at Q'^/M'^ — oo in a way which guarantees smoothness [14] [15], but it seems 
phenomenologically motivated, with no strict definition of the ordering and certainly no expressions 
for parton distributions and coefficient functions in the intermediate region. The ACOT group have 
a prescription which involves switching from n/ to n/ + 1 massless flavours in the evolution, and 
a way of determining the VFNS coefficient functions [22] [23] which at low orders appears to be 
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the same as prescribed in (3.5)-(3.9). However, their way of eUminating the free choices in the 
heavy quark coefficient functions involves assuming that the behaviour is as if there is intrinsic 
charm in the proton at all scales above the transition point, rather than charm being generated 
almost entirely from the gluon. This leads to coefficient functions having thresholds in = 
rather than = 4M^ , and a mixing of orders being required (and a complicated renormalization 
scale being advantageous) in order to ensure cancellations and that smooth behaviour occurs, e.g. 
the 0{as{n'^)) gluon coefficient function must appear at the same time as the zeroth order quark 
coefficient function. This mixing of orders is incorrect, being at odds with well-ordered asymptotic 
expressions, but removing it results in a lack of smoothness in the structure functions. Even when 
this mixing is retained the behaviour of the longitudinal structure function is still not smooth. 
The MRRS procedure [28] is based on the leading log limit of Feynman diagrams rather than the 
renormalization group and as such incorporates mass-dependent effects in the evolution, but seems 
more difficult to define formally to all orders in ccg. The definition of the heavy quark coefficient 
functions uses similar reasoning to ACOT, but in this case with ordering such that it reduces to 
the correct well-ordered form asymptotically. These coefficient functions along with the imposition 
of this correct ordering leads to an unphysical lack of smoothness in the structure functions (which 
is made slightly worse by the mass-dependent contributions to the evolution), particularly for the 
longitudinal structure function. Our prescription has none of the above problems. It is well-defined 
to all orders, reduces to correct well-ordered expressions at both low and high Q^, and exhibits 
precisely the behaviour one would expect. Hence, we believe that our prescription is the best 
currently available to describe the heavy quark contribution to structure functions. 

Before finishing let us mention a couple of points in which our treatment is incomplete. Firstly, 
we have assumed that there is no intrinsic charm in the nucleon. Eqn. (2.1) is formally correct up to 
the quoted error, but this error has an unknown numerical factor and may be enhanced by functions 
of X. It appears that for intrinsic charm the numerical factor of this "higher twist" correction is 
rather large and that the contribution is enhanced by a factor of (1 — a;)~^. Therefore, at large x , 
where the leading twist contribution to the charm structure function is not large anyway, it seems 
as though the "higher twist" intrinsic charm may constitute an important part of the total charm 
structure function [35]. The treatment of this correction is outside the scope of this paper, which 
deals only with the "leading twist" contribution to the structure function, and we believe it is not 
naturally dealt with in any other VFNS. However, it seems very unlikely that in most of the region 
where there is current data on the charm structure function, or where the charm contribution is 
a sizable fraction of the total structure function, that this "higher twist" contribution plays any 
significant role at all. Using the type of values expected for this intrinsic charm (see e.g. [36]) then 
adding to our values does bring the x = 0.422 prediction in line with the EMC data point, raises 
the X = 0.237 predictions quite significantly (but neither really helps or hinders the comparison to 
the three data points), raises the x = 0.133 predictions a little (tending to make the comparison a 

40 



little worse), and has negligible effect for lower x. Hence, the x = 0.422, = 78GeV^ EMC data 
point may be seen as some evidence for this "higher twist" intrinsic charm. 

Finally we note that throughout this paper we have completely ignored the problem of en- 
hancement of higher orders in by ln(l/x) terms. These terms certainly do have the potential to 
quantitatively alter the results of this paper. Correctly including the leading ln(l/x) terms within 
the context of only massless quarks is a complicated procedure, though it does appear to improve 
the description of small x data [20] . Some results on heavy quark coefficient functions which include 
leading ln(l/x) terms already exist [37][38]. It would clearly be desirable to extend this work and 
to include both the correct treatment of leading ln(l/x) terms and a correct description of heavy 
quark results within a single framework. Work along these lines is currently in progress. 

Appendix A. Charged Current Structure Functions. 

The treatment of the charged current structure function follows exactly the same reasoning as 
for the neutral current case. Let us consider F2{x, Q^). Equations (3.5)-(3.9) are derived in exactly 
the same way, but now take a different form because there are no nonsinglet coefficient functions. 
For the case where a heavy quark is produced directly by the interaction with the W boson, which 
we call the heavy quark structure function, we have 

^Hg - ^gg,H ® <-^Hg + ^f\g,H ® '-'Hq + ^Hg ® ^HH ' {^■^) 

and 

C'h^'' = ^ C]^r' + Kq% + nfA^,%] cYg^' + 4,,, ^ cVg^'. (A.2) 

We note that what we have denoted the charm quark structure function here may be interpreted 
physically as the unlike sign dimuon contribution. In the case where the W boson directly produces 
a light quark, which we call the light quark structure function, we have 

Cig + C^f/'" = ^ Cl^^' + nMf,,^ ^ C^J-^^ + ^ C,^/'^^ (A.3) 

and 

riPS r'FF,PS ^ iPS ^r'VF,PS , aPS ^ r^VF,PS . aS ^ riVF,s /a a\ 

It is not only the absence of the nonsinglet coefficient functions which is different, the ordering of 
the other coefficient functions also changes, in particular C2 g.g. , i ^ j, begins at zeroth order. This 
changes the form of the relationship between the FFNS and the VFNS coefficient functions. For 
example, examination of (A.l) reveals that we have the trivial equality 

C^^jiz, Q^M^) ^ C^^jiz, Q'/M^ (A.5) 
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whereas now we have the nontrivial relationship 



QVM^) = Cl^^iz, Q'/M') - iHQ'm + c..)P°, (A.6) 

e.g. the zeroth order coefficient function for a charm quark to interact with a W~ to produce a 
strange quark is undetermined. As in the previous case we determine this zeroth order heavy quark 
coefficient function by demanding continuity of the In Q^-derivativc of tlic structure function, in 
the gluon sector (again at lowest order we have complete continuity), along with demanding the 
correct asymptotic result. Unlike the neutral current case this time it is the strange quark (or 
down quark) structure function on which the condition is imposed, rather than the charm quark 
structure function. This is because at lowest order the charm quark structure function is completely 
independent of the charm quark distribution, whereas the light quark structure functions do depend 
on it. However, in complete analogy with the neutral current case our constraint results in 

Ci:r^QVM^) ^ = '""'''^{^^^^^^^^^ (A.7) 

where the left hand side automatically has the correct threshold behaviour and the right-hand side 
Pgg{z) as Q^/M^ DO. Using this explicitly in (A.6) then results in the C^^'^(2:, QV-^^) 
reducing to the correct massless MS limit as Q^/M"^ — ^ oo, as it must by construction. 

This procedure can be repeated at all orders in exactly the same way as for the neutral 
current structure function. This time there are only two coefficient functions to be determined, 
C^Hq^i^i Q"^ /M'^) as we have just seen, and which exists at all orders, and C^H^^iz, Q^/M^). The 
later begins at 0{al{Q'^)) and will be determined by demanding continuity of the In Q^-derivative 
of the structure function where a heavy quark is produced directly at the interaction vertex of the 
W boson at 0{al{Q'^)). The extension to the longitudinal charged current structure function is 
also easily achieved using the above results and the discussion of the longitudinal neutral current 
structure function in S7. 
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Table 1. 



Prescription for the order by order in asiQ"^) determination of the VFNS coefficient functions 
for F2(x,Q2). 



Order of Eqn. 
equality 

aO(Q2) (3.7 

3.5 
3.7 
3.8 

3.5 
3.6 
3.7 
3.8 
3.9 

3.5 
3.6 
3.7 
3.8 
3.9 

3.5 
3.6 
3.7 
3.8 
3.9 



Coefficient functions determined 



iVF.NS.O 



CIhS"" (by continuity of (S#)m^ at 0(a,(Q2))), ^--^^ 



^VF,l 



f-iVF,NS,l 

r<VF,i 

^VF,S,1 



CIhh^ (by continuity of (^§^)m2 in gluon sector at 0{al{Q^))), C^^Hg 

C, 

a 



VF,2 



VF,2 
2,Hq 
VF,NS,2 
2,qq 



a 



'2,qg 
VF,PS,2 



2,qq 

^VF.S,2 



CIhh^ (by continuity of (^^^)m2 in gluon sector at C>(a^(Q^))), C^^'g 

a 
a 



VF,3 



VF.3 
2,Hq 
VF,3 

iVF,2 



ClqH (by continuity of (^^|^)m2 in gluon sector at 0{al{Q'^))), 



VF,3 
19 



r<VF,PS,3 
^2,qq 



(by continuity of {^)m- in gluon sector at 0«(g2))), Cl^j^^ 



dF2,H 



VF,n 



riVF,n 
^2,Hq 
riVF,n 
'^2,qq 
iVF,n-l 



Cl^R-' (by continuity of (^)m^ in gluon sector at 0{a^{Q^))), Cl^^ 



VF,n 



r<VF,PS,n 
^2,qq 
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Table 2. 



Comparison of quality of fits for a wide variety of structure function data [2] [3] [33] [39] [40] [41] using 
our prescription for heavy quarks at NLO (TR) and the NLO fits MRRS and MRS(R2). We do not 
include the small-x, low-Q^ data in the second of [2] in our fit in order to make a direct comparison 
with the previous fits. 



Experiment 


data 










points 


TR 


MRRS 


MRS(R2) 


HI F^" 


193 


135 


133 


149 


ZEUS 


204 


274 


290 


308 


BCDMS F^P 


174 


262 


271 


320 


NMC F^P 


130 


144 


145 


135 


NMC F^'^ 


130 


112 


119 


99 


E665 F^^ 


53 


61 


60 


62 


E665 F^"^ 


53 


53 


54 


60 


SLAC F^^ 


70 


98 


96 


95 
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Table 3. 



Prescription for the order by order in asiQ"^) determination of the VFNS coefficient functions 
for FL(x,g2). 



Eqn. 



Coefficient functions determined 



3.5 
3.7 
3.8 

3.5 
3.6 
3.7 
3.8 
3.9 

3.5 
3.6 
3.7 
3.8 
3.9 

3.5 
3.6 
3.7 
3.8 
3.9 



riVF,i 

iVF,S,l 



CL^hh" (by continuity of {§^)m^ in gluon sector at Cl^j}; 



VF,2 



VF,NS,2 



^L,qq 
^VF,2 



c 



'L,qg 
VF,PS,2 



L,qq 
iVF,S,2 



Cl'ii'ii' (by continuity of (^)m^ m gluon sector at 0{ai{Q''))), Cl^^^ 



VF,3 



riVF,z 

^L,Hq 
^VF,3 
^L,qq 
iVF,2 



^UqH (by continuity of {j^^)m^ in gluon sector at 0{al{Q'^))), 



VF,3 
19 



riVF,PS,z 

'^L,qq 



Cl%^H~^ (by continuity of (|^)m^ in gluon sector at 0«(g2))), d'^^^^ 



■ dFL,H 



VF,n 



/^VF,n 
L,Hq 

^L,qq 
^VF,n-l 



Cl!qH~' (by continuity of (^^)m^ in gluon sector at 0{a^,{Q^))), 



C 



VF,PS,n 
L,qq 



VF,n 
19 



In each case C^^'^(z, M"^ /Q^) is determined by introducing the asymptotic form multiplied 
by (1 — M'^/Q'^)v and multiplying the terms determined by continuity by M'^/Q'^. 
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Figure Captions 



Fig. 1. Charm quark structure function, F2.c{x,Q'^) for x = 0.05 and x = 0.005 calculated using 
the ACOT "LO" prescription, our input parton distributions evolved at LO and renor- 
malization scale /i^ = Q^. Shown are the total, the two contributions due to convolution 
of the coefficient function C^j^'g (z, Q^/m^) with the gluon distribution (the subtraction 
term making a negative contribution), and the contribution directly due to the charm 
quark. 

Fig. 2. Same as fig. 1, but with the partons evolved at NLO. 

Fig. 3. Charm quark structure function, F2^cix,Q'^) for x = 0.05 and x = 0.005 calculated using 
our LO prescription, our input parton distributions evolved at LO and renormalization 
scale /U^ = . Also shown are the continuation of the LO FFNS expression and the ZM- 
VFNS expression both calculated using the same parton distributions and same choice of 
scale. 

Fig. 4. Same as fig. 3, but with NLO prescriptions and NLO parton distributions. 

Fig. 5. Our prediction for F2^c(,x,Q'^) using our NLO prescription, the NLO partons obtained 
from our global fit and three different values of rric compared with the EMC and HERA 
data. 

Fig. 6. The ratios -F2,c/^2 and .^2,6/^2 at fixed values of resulting from our NLO parton 
distributions and taking rric = 1.35GeV and mi, = 4.3GeV. The experimental data point 
shows the estimate from ref. [4] for -F2,c/-P2 in the kinematic range lOGeV^ < Q2 < 
lOOGeV^ 

Fig. 7. Charm quark structure function, Fl^c{x,Q'^) for x = 0.05 and x = 0.005 calculated using 
our LO prescription, our input parton distributions evolved at LO and renormalization 
scale fP' = Q^. Also shown are the continuation of the LO FFNS expression and the ZM- 
VFNS expression both calculated using the same parton distributions and same choice of 
scale. 

Fig. 8. Same as fig. 7, but with NLO prescriptions and NLO parton distributions. 

Fig. 9. Our prediction for i?(x, Q^) using our NLO prescription, the NLO partons obtained from 
our global fit and rric = 1.35GeV compared with the NMC data [33]. Also shown is 
the prediction obtained using the same parton distributions but for the NLO ZM-VFNS 
prescription. The curves are computed using g2 = L3GeV^ for x < 0.0077 and = 
262x109 for x > 0.0077. 
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